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1 Introduction 

In these lecture notes, corresponding roughly to lectures given at the summer school in Chia La- 
guna, Italy, in September 1997, an attempt is made to present the physics of three-dimensional 
interacting fermion systems (very roughly) and that of their one-dimensional counterparts, the 
so-called Luttinger liquids (in some more detail). These subjects play a crucial role in a num- 
ber of currently highly active areas of research: high temperature and organic superconductors, 
quantum phase transitions, correlated fermion systems, quantum wires, the quantum Hall effect, 
low-dimensional magnetism, and probably some others. Some understanding of this physics thus 
certainly should be useful in a variety of areas, and it is hoped that these notes will be helpful in 
this. 

As the subject of these lectures was quite similar to those delivered at Les Houches, some 
overlap in the notes|jl[ was unavoidable. However, a number of improvements have been made, 
for example a discussion of the "Klein factors" occurring in the bosonization of one-dimensional 
fermions, and new material added, mainly concerning spin chains and coupled Luttinger liquids. 
Some attempt has been made to keep references up to date, but this certainly has not always been 
successful, so we apologize in advance for any omissions (but then, these are lecture notes, not a 
review article). 



2 Fermi Liquids 

Landau's Fermi liquid theory[0, ^ Q is concerned with the properties of a many-fermion system 
at low temperatures (much lower than the Fermi energy) in the normal state, i.e. in the absence 
or at least at temperatures above any symmetry breaking phase transition (superconducting, 
magnetic, or otherwise). The ideal example for Landau's theory is liquid helium 3, above its 
superfiuid phase transition, however, the conceptual basis of Landau's theory is equally applicable 
to a variety of other systems, in particular electrons in metals. Quantitative applications are 
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however more difficult because of a variety of complications which appear in real systems, in 
particular the absence of translational and rotational invariance and the presence of electron- 
phonon interactions, which are not directly taken into account in Landau's theory. Subsequently, 
I will ffist briefly discuss the case of a noninteracting many-fermion system (the Fermi gas), and 
then turn to Landau's theory of the interacting case (the liquid), first from a phenomenological 
point of view, and then microscopically. A much more detailed and complete exposition of these 
subjects can be found in the literature p, ^ ^ |], ||. 

2.1 The Fermi Gas 

In a noninteracting translationally invariant systems, the single-particle eigenstates are plane 

waves 

Ik) ^ (2.1) 
with energy 

= (2.2) 

where fl is the volume of the system, and we will always use units so that h = 1. The ground 
state of an A^-particle system is the well-known Fermi sea: all states up to the Fermi wavevector 
kp are ffiled, all the other states are empty. For spin-1/2 fermions the relation between particle 
number and kp is 

N = n'^ . (2.3) 

The energy of the last occupied state is the so-called Fermi energy E-p = kp^ / (2m), and one easily 
verifies that 

^,.f!^.,(r^O) (2.4) 

i.e. E-p is the zero-temperature limit of the chemical potential {Eq{N) in the formula above is the 
ground state energy). 

It is usually convenient to define the Hamiltonian in a way so that the absolute ground state has 
a well-defined fixed particle number. This is achieved simply by including the chemical potential 
/i in the definition of the Hamiltonian, i.e. by writing 

i/ = E^knk , (2.5) 

k 

where nk is the usual number operator, ^k = £^k — /W, and the spin summation is not written 
explicitly (at finite temperature this of course brings one to the usual grand canonical description 
where small fluctuations of the particle number occur). With this definition of the Hamiltonian, 
the elementary excitations of the Fermi gas are 

• addition of a particle at wavevector k ((5nk = !)• This requires |k| > kp, and thus the energy 
of this excitation is ^k = Ck — > 0. 

• destruction of a particle at wavevector k (5rak = —1), i-e. creation of a hole. This requires 
\k\ < kp, and thus the energy is 6:^ = /i — ek > 0. 
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(a) 



(b) 



Figure 1: (a) the energy-momentum relation for the elementary particle {k > kp) and hole {k < k-p) 
excitations; and (b) the particle-hole continuum. 

The dispersion relation of the elementary particle and hole excitation is shown in fig.|l]a. These 
excitations change the total number of particles. Construction of states at constant particle 
number is of course straightforward: one takes one particle from some state k, with \k\ < k-p, and 
puts it into a state k', with |k'| > kp. These particle-hole excitations are parameterized by the 
two quantum numbers k,k' and thus form a continuum, as shown in fig.^. The restriction on the 
allowed values of k, k' insures that all particle-hole states have positive energy. Higher excited 
states, i.e. states with many particles and many holes, are straightforwardly constructed, the only 
restriction being imposed by the Pauli principle. 

Thermodynamic quantities are easily obtained and are all determined by the density of states 
at the Fermi energy. For example, the specific heat obeys the well known linear law C{T) = 7T, 
with 

7 = ^NiEp)kl (2.6) 

and similarly the (Pauli) spin susceptibility x ci^nd the compressibility k are given by 

X = 2N{Ep)fil (2.7) 
K = 2N{Ep)lp^ . (2.8) 



Here for the quadratic dispersion relation (|2.2|) the density of states (per spin) at the Fermi energy 



is given by N{Ep) = mkp / {2tt'^) , but it should be emphasized that eqs. (|2.6| ) to ( p.8| ) are valid for 



an arbitrary density of states, in particular in solids where bandstructure effects can change the 
electronic dispersion relation quite drastically. Thus, for noninteracting electrons one expects the 
so-called "Wilson ratio" 



Rw = (2.9) 
3/iB 7 



to be unity, independently of details of the bandstructure. Any deviation from unity is necessarily 
an indication of some form of interaction effect. 
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2.2 Landau's theory of Fermi Liquids 
2.2.1 Basic hypothesis 

Landau's theory is to a large extent based on the idea of a continuous and one-to-one correspon- 
dence between the eigenstates (ground state and excited states) of the noninteracting and the 
interacting system. For this to be an acceptable hypothesis it is crucial that the interactions do 
not lead to any form of phase transition or symmetry-broken ground state. 

In particular one can consider a state obtained by adding a particle (with momentum |p| > /cf) 
to the noninteracting ground state: 



state of the noninteracting system. Now we add some form of particle-particle interaction. In 
a translationally invariant system, interactions conserve total momentum, and thus even after 
switching on the interaction the state still has total momentum p. However, the interaction of 
the added particle with the filled Fermi sea, as well as the interaction of the particles in the 
sea amongst themselves, will change the distribution of particles in k-space, and of course also 
modify the energy of our state. The complex formed by the particle added at p and the perturbed 
distribution of the other particles is called a Landau quasiparticle. The Pauli principle implied 
|p| > kp in the absence of interactions, and by the continuity hypothesis the same restriction 
remains valid in the interacting case. In particular, the value of kp, which imposes a lower limit 
on the allowed momentum of the quasiparticle, is unchanged by the interactions. 

Analogous considerations can be performed for a state obtained by destruction of a particle 
(e.g. creation of a hole): 



Note that due to the momentum — p the total momentum of this state is indeed p. 

The quasi-particle concept has a certain number of limitations, mainly due to the fact that, 
as will be discussed below, the lifetime of a quasi-particle is finite. However, for excitations close 
to kp one has l/r oc (e — Ep)^, i.e. the lifetime becomes much longer than the inverse excitation 
energy, and the quasi-particles therefore are reasonably well defined. In practice, this means that 
Landau's theory is useful for phenomena at energy scales much smaller than the Fermi energy, but 
inapplicable otherwise. In metals, where E'p ~ 3 . . . 5eV, this restriction is not too serious when 
one is concerned with thermodynamic or transport properties. One should also note that the 
ground state energy itself has important contributions from states well below Ep, and therefore is 
not accessible to Landau's theory. 

2.2.2 Equilibrium properties 

In order to derive physical quantities from the picture of the low-energy excitations, we need some 
information about the energetics of the quasiparticles and of their interactions. To be specific, 
starting from the ground state quasiparticle distribution 

no(k) = 1 if |k| < kp 





(2.11) 



= if Ikl > kp 



(2.12) 
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one considers changes in quasiparticle occupation number of the the form nQ{k) nQ{k) +5n{k), 
i.e. Sn{k) = 1 represents an excited quasi-particle, 6n{k) = —1 an excited quasi-hole (with the 
notation k = (k, a), and o" =1, | the spin index). The corresponding change in energy is 

6E = Y.elSnik) + ^J2fik^k')Hk)6n{k') , (2.13) 
k kk' 

where the first and second term represent the energy of a single quasi-particle and the interaction 
between quasiparticles, respectively. To be more precise, we assume that the chemical potential is 
included in the Hamiltonian, as in eq.( p.5|) . Consequently, vanishes on the Fermi surface, and, 
given that we are mainly interested in phenomena in the vicinity of k-p, it is sufficient to retain 
the lowest order term in an expansion around |k| = k-p. One thus writes 

^k=%(|k|-M , (2.14) 

m* 

thus defining the effective mass m* which is different from the "bare" mass m due to interaction 
effects that could in principle be calculated from a microscopic theory of the system. 

The energy of a quasi-particle added to the system is easily obtained from eq.( p.l3|) by cal- 
culating the difference in SE between a state with 6n{k) = 1 and a state with 6n{k) = 0. One 
finds 

ek = ei + Wf{k,k')6n{k') , (2.15) 
k' 

i.e. the energy of an added quasi-particle is not just the "bare" quasiparticle energy e:^ but also 
depends, via the interaction term, on the presence of the other quasi-particles. Given that the non- 
interacting particles obey Fermi-Dirac statistics, the quasi-particles do so too, and consequently, 
the occupation probability of a quasi-particle state is given by 

nik) = . (2.16) 

Note that the full and not the bare quasi-particle energy enters this expression. In principle, n{k) 
thus has to be determined self-consistently from eqs. (|2.15| ) and (|2.16|) . 

For the subsequent calculations, it is convenient to transform the quasiparticle interaction 
f{k,k'). First, spin symmetric and antisymmetric /-functions are defined via 

/(kt,k't) = r(k,k') + r(k,k') 

/(kT,k'i) = r(k,k')-r(k,k') (2.17) 

Moreover, given the implicit restrictions of the theory, one is only interested in processes where all 
involved particles are very close to the Fermi surface. Under the assumption that the interaction 
functions are slowly varying as a function of k, one then can set |k| = |k'| = kp. Because of 
rotational symmetry, the /-functions then can only depend on the angle between k and k', called 
6. One can then expand the /-function in a Legendre series as 

oo k ■ k' 

/'^'^(k,k') = Y: frPdcosO) ,cose = , (2.18) 

L=0 '^F 

where the Pl are the Legendre polynomials. Finally, one usually puts these coefficients into 
dimensionless form by introducing 

Fl'' = '^ir ■ (2.19) 
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We are now in a position to calculate some equilibrium properties. The first one will be the 
specific heat at constant volume 



Cc. 



(2.20) 



where U is the internal energy. The temperature-dependent part of U comes from thermally 
excited quasi-particles, as determined by the distribution ( p.l6| ). In principle, in this expression 
Ek is itself temperature-dependent, because of the temperature dependent second term in eq.( |2.15D . 
However, one can easily see that this term only gives contributions of order T^, and therefore can 
be neglected in the low-temperature limit. Consequently, one can indeed replace Ek by e^, and 
then one only has to replace the bare mass by m* in the result for a non-interacting system to 
obtain 



Cc 



(2.21) 



The spin susceptibility (at T = 0) is related to the second derivative of the ground state energy 
with respect to the (spin) magnetization M: 



X 



n 



f2.22) 



Spin magnetization is created by increasing the number of | spin particles and decreasing the 
number of | spins (M = fiBi^i ~ ^i))) i-^- t>y changing the Fermi wavevectors for up and down 
spins: kp ^ kp + 6kp for a =t and kp ^ kp — 5kp for a =[. 

By calculating with eq.( p.l3| ) the corresponding change of the ground state energy, we obtain 
from eq.(|2]2|): 



X 



fi'^kpm* 



1 + 7r2 



(2.23) 



Note that here, and contrary to the specific heat, interactions enter not only via m* but also ex- 
plicitly via the coefficient Fq, which is the only coefficient that appears here because the distortion 
of the Fermi distribution function is antisymmetric in the spin index and has spherical symmetry 
(L = 0). The Wilson ratio is then 



R 



1 



w 



1 + F,^ 

Following a similar reasoning, one can calculate the compressibility k of a Fermi liquid: 



(2.24) 



K = — 



1 dn 



m*kp 



7rV(l + ^0 



(2.25) 



It is also interesting that in a translationally invariant system as we have considered here, 
the effective mass is not independent of the interaction coefficients. One can show indeed, by 
exploiting the Galilean invariance of the system, that 



— = l + Fl/3 . 
m 



(2.26) 
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2.2.3 Nonequilibrium properties 

As far as equilibrium properties are concerned, Landau's theory is phenomenological and makes 
some important qualitative predictions, the most prominent being that even in the presence of 
interactions the low-temperature specific heat remains linear in temperature and that the spin 
susceptibility tends to a constant as T — > 0. However, Landau's theory has little quantitative 
predictive power because the crucial Landau parameters have actually to be determined from 
experiment. The situation is different for non-equilibrium situations, where the existence of 
new phenomena, in particular collective modes, is predicted. These modes are another kind of 
elementary excitations which, contrary to quasiparticles, involve a coherent motion of the whole 
system. We shall not enter into the details of the treatment of non-equilibrium properties of Fermi 
liquids (see refs.[^, ||, ||) and just briefly sketch the general conceptual framework and some of the 
more important results. To describe non-equilibrium situations, one makes two basic assumptions: 

• Deviations from equilibrium are described by a Boltzmann equation for a space- and time- 
dependent quasiparticle distribution function n(k, r, t), which describes the density of quasi- 
particles of momentum and spin (k, a) at point r and time t. At equilibrium, n is of course 
given by eq. (|2.12| ). The fact that in the distribution function one specifies simultaneously mo- 



mentum and position of course imposes certain restrictions, due to the quantum-mechanical 
nature of the underlying problem. More precisely, spatial and temporal variations of the 
distribution function have to be slow compared to the typical wavelength and frequency of 
the quasiparticles. We have then the conditions frlql, < -^f, where q and |a;| set the 
scale of the spatial and temporal variations of n{r,t). 



• Because of the r-dependent n, the quasiparticle energy is itself, via eq. (|2.15|) , r-dependent. 
One then assumes the following quasi-classical equations of motion 

r = Vk£:k(r) 

k = -Vr£k(r) . (2.27) 

Note in particular that a space-dependent distribution function gives rise, via the f{k, k') 
function, to a force acting on a quasiparticle. 

By linearizing the Boltzmann equation and studying the coUisionless regime, where the collision 
term in the Boltzmann equation can be neglected, one finds collective mode solutions which 
correspond to oscillations of the Fermi surface. The most important one is the longitudinal 
symmetric mode which, like ordinary sound, involves fluctuations of the particle density. This 
kind of sound appears, however, in a regime where ordinary sound cannot exist (the existence of 
collisions is indeed crucial for the propagation of ordinary sound waves) and is a purely quantum 
effect. Since collisions can always be neglected at very low temperatures, this new kind of sound 
has been called by Landau zero sound. The collision term of the Boltzmann equation is on the 
contrary essential to calculate the quasiparticle lifetime r. One can find indeed, for a quasiparticle 
of energy e-p 

r-^ oc ^*3 (^^)' + 4 _ (2.28) 

The most important result here is the divergence of the lifetime for low energies and temperatures 
as r oc max(£:p, T)~^, so that the product SpT in fact diverges as the Fermi surface is approached. 
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This shows that the quasiparticle becomes a well-defined (nearly-) eigenstate at low excitation 
energies, i.e. in the region where Landau's theory is applicable. On the other hand, at higher 
energies the quasiparticle becomes less and less well-defined. One may note that initially we had 
assumed that a quasiparticle is an exact eigenstate of the interacting system, which was obtained 
from a noninteracting eigenstate by switching on the interaction, and therefore should have infinite 
lifetime. We now arrive at the conclusion that the lifetime is not strictly infinite, but only very 
long at low energies. In the following section we will try to clarify this from a microscopical point 
of view. 

2.3 Microscopic basis of Landau's theory 

At our current knowledge, it does not seem generally possible to derive Landau's theory starting 
from some microscopic Hamiltonian, apart possibly in perturbation theory for small interactions. 
It is however possible to formulate the basic hypotheses in terms of microscopic quantities, in 
particular one- and two-particle Green functions. This will be outlined below. 

2.3.1 Quasiparticles 

As far as single particle properties are concerned it is sufficient to consider the one-particle Green 
function 



where r is the usual (Matsubara) imaginary time. In this quantity, interaction effects appear via 
self-energy corrections S in the Fourier transformed function 



Here is the bare particle energy, without any effective mass effects. Excitation energies of the 
system then are given by the poles of Q{k,uj). In these terms. Landau's assumption about the 
existence of quasiparticles is equivalent to assuming that u) is sufficiently regular close to the 
Fermi surface as to allow an expansion for small parameters. Regularity in (k, u) space implies 
that in real space the self-energy has no contributions that decay slowly in time and/or space. 
Given that the self-energy can be calculated in terms of the effective interaction between particles 
this is certainly a reasonable assumption when the particle-particle interaction is short-range 
(though there is no formal prove of this). For Coulomb interactions, screening has to be invoked 
to make the effective interaction short ranged. 

One can further notice that ll{kp,0) just renormalizes the chemical potential. Given that we 
want to work at fixed particle number we can absorb this term in the effective fi. Expanding then 
to first order around the Fermi surface, the Green function takes the form 




(2.29) 




(2.30) 




(2.31) 



where £^ has the form ( p.l4|) of the phenomenological approach, with 
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Figure 2: The momentum distribution function 
acting (dashed line) cases. 



(aj^Ok) in the interacting (full line) and noninter- 



and the quasipartide renormalization factor is 

-1 



z 



1 - 



duj 



All derivatives are to be taken at the Fermi surface and at uj 
rule imposes that the frequency-integrated spectral density 



(2.33) 



0. One should notice that a sum 



A{k,uo) 



TT 



-\m.Q{k,iuj uj + i6) 



(2.34) 



equals unity. Consequently, in order to fulfill the sum rule, if z < 1 there has to be a contribution 
in addition to the quasipartide pole in eq.( p.31|) . This is the so-called "incoherent background" 
from single and multiple particle-hole pair excitations which can extend to rather high energies 
but becomes small close to the Fermi surface. 

The form ( p.31| ) gives rise to a jump in the momentum distribution function at kp of height 
z, instead of unity in the noninteracting case (fig.0). In addition to the jump, the incoherent 
background gives rise to a contribution which is continuous through kp. 

A finite quasipartide lifetime arises if the expansion of S(fc, u) is carried to second order. Then 
eq. (p.31|) generalizes to 



Q{k,uj) 



iuj — e^^ + isign(co')r(co') ^ 



(2.35) 



where t{uj) is typically given by an expression like eq. (|2.28|) . 



2.3.2 Quasipartide interaction 

The quasipartide interaction parameters f{k,k') are expected to be connected to the the two- 
particle vertex function. This function, which we will denote T^'^\Pi, P2; K) describes the scatter- 
ing of two particles from initial state Pi, P2 to the final state Pi — i^, P2 + K, and the notation is 
Pi = {pi,uJi, ai). The contribution of first and second order in the interaction potential V{\i) are 
shown in fig.^. Let us now study the case of small transfer /T, but arbitrary Pi 2, only restricted to 
be close to the Fermi surface. One then notices that diagram |^ (part c) gives rise to singularities. 
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because the poles of the two intervening Green functions coalesce. On the other hand diagrams 
^ (part a) and ^ (part b) remain nonsingular for small K. This motivates one to introduce a 
two-particle irreducible function F*^^) which is the sum of all contributions which do not contain a 
single product Q{Q)Q{K + Q). This function then is nonsingular for small K, and consequently 
the total vertex function is determined by the integral equation 



T{P,,P^-K) = f{P,,P^) 



dijd'qt{Pi, Q)g{Q)g{K + g)r(Q, P2, k) 



(2.36) 



For simplicity, the spin summation is omitted here. The singular contribution now comes from 
small K and Q in the vicinity of the Fermi surface. In this area the Q-dependence of the F's in 
eq.( p.36|) is non-singular and can be neglected. The energy and radial momentum integral over Q 
can then be done, leading to 



F(Pi,P2;i^) = f(Pi,P2) + 



(27r)3 



d''n^f{p,,Q)- 



q ■ k 



tipq ■ k 



nQ,P2,K) , 



(2.37) 



where q is a vector on the Fermi surface, and d^Qq is the corresponding angular integration. Here 
only the quasiparticle pole in Q has been taken into account. The contribution from the incoherent 
parts can in principle be absorbed into the definition of F. 

The expression (|2.37| ) is clearly singular because it has radically different behavior according 
to whether one first sends k or a; to zero. In the first case, the limit of F can be related to 
the Landau /-function, while the second one is relevant for the calculation of the transition 
probabilities determining the lifetime of quasiparticles. Here we will consider only the former 
case. Sending k to zero first one finds straightforwardly 



limflim F) 

a;-»0 k^O 



F-(Pl,P2)=f(Pi,P2) 



(2.38) 



Closer inspection then shows that in this case the poles of the two Green functions in eq. (|2.36| ) are 
always on the same side of the real axis and consequently the singular term in eq.( 2.36|) vanishes. 
To make the identification between F'^(Pi,P2) and the Landau /-function we notice that the 
density response function at energy-momentum K, whose poles give the collective (zero-sound) 
modes, contains the interactions via F(Pi, P2, i^'). In particular, the existence of a pole in the 
response function implies a pole in F. The comparison between the equations for this pole and 



Pi 


Pi5 




(2) 

r (PpPjK) 




> 







P1+P2-Q 



V(k) + 




Po+K 





V(k-p, +pj 



Q 
(a) 

PpP^K+Q 




Q 

(b) 



Qi I i^+Q 




(a) 



(b) 



Figure 3: The lowest (a) and second order (b) contributions to the two-particle vertex F^^^. Note that 
the external lines do not represent actual Green functions but only indicate the external "connections" 
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those which can be obtained through the Boltzmann equation within Landau's theory allows then 
the identification 



f{k,k') = z^T'^{k,k') . (2.39) 
2.4 Summary 

The basic assumption of Landau's theory is the existence of low-energy quasiparticles with a very 
long lifetime, and their description in terms of a rather simple energy functional, eq. (|2.13|) . From 
this a number of results for thermodynamic properties is obtained. At this level, the theory is 
of little quantitative power because the Landau parameters are not determined. Qualitatively, 
however, the predictions are important: the low-temperature thermodynamic properties of an 
interacting fermion system are very similar to those of a noninteracting system, the interactions 
only lead to quantitative renormalizations. Actual quantitative predictions are obtained when 
one extends the theory to nonequilibrium properties, using the Boltzmann equation.]^ A new 
phenomenon predicted (and actually observed in ^He |]lO|) is the existence of collective excitations, 
called "zero sound" . This approach also allows the calculation of the quasiparticle lifetime and its 
divergence as the Fermi energy is approached, as well as the treatment of a number of transport 
phenomena. 

As already mentioned, the ideal system for the application of Landau's theory is ^He, which 
has both short-range interaction and is isotropic. The application to electrons in metals is more 
problematic. First, the interactions are long-ranged (Coulombic). This can however be accom- 
modated by properly including screening effects. More difficulties, at least at the quantitative 
level, arise because metals are naturally anisotropic. This problem is not of fundamental nature: 
even when the Fermi surface is highly anisotropic, an expansion like eq.( p.l3D can still be written 



down and thus interaction parameters can be defined. However, a simple Legendre expansion like 
eq.( p.l8|) is not in general possible and the description of the quasiparticle interaction in terms of a 



few parameters becomes impossible. An exception case, with a very nearly spherical Fermi surface, 
are the alkali metals, where a determination of Landau parameters can indeed be attempted. It 
should be noticed that the difficulties with the Landau description of metals are not of conceptual 
nature and in particular do not invalidate the quasiparticle concept but are rather limitations on 
the usefulness of the theory for quantitative purposes. 

Landau's theory can be interpreted in terms of microscopic quantities like Green functions (the 
quasiparticle pole) and interaction vertices, as discussed above. It should however be emphasized 
that these arguments do provide a microscopic interpretation of Landau's picture, rather than 
proving its correctness. Similar remarks apply to the calculated diverging quasiparticle lifetime: 
this at best show that Landau's picture is internally consistent. Considerable progress towards a 



deeper formal understanding of Fermi liquid theory has been made in recent years. 1 11, 12 



3 Renormalization group for interacting fermions 

In this chapter, we will consider properties of interacting fermions in the framework of renor- 
malization group theory. This will serve two purposes: first, the treatment of one-dimensional 
interacting fermions, which will be considered in considerable detail in the following chapters, gives 
rise to divergences which can only be handled by this approach. Results obtained in this way will 
be an essential ingredient in the subsequent discussion of "Luttinger liquids" . More generally, the 
renormalization group method will clarify the status of both Landau's Fermi liquid theory and 
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the Luttinger liquid picture as renormalization group fixed points, thus estabhshing a hnk with 
a number of other phenomena in condensed matter physics. We will formulate the problem in 
terms of fermion functional integrals, as done by Bourbonnais in the one-dimensional case 



and more recently for two and three dimensions by Shankar ^M. For the most part, I will closely 
follow Shankar's notation. 

Before considering the interacting fermion problem in detail, let us briefiy recall the general 
idea behind the renormalization group, as formulated by Kadanoff and Wilson: one is interested 
in the statistical mechanics of a system described by some Hamiltonian H. Equilibrium properties 
then are determined by the partition function 

Z= E e~^" = E , (3.1) 

configurations configurations 

where the second equality defines the action S = [3H. Typically, the action contains degrees of 
freedom at wavevectors up to some cutoff A., which is of the order of the dimensions of the Brillouin 
zone. One wishes to obtain an "effective action" containing only the physically most interesting 
degrees of freedom. In standard phase transition problems this is the vicinity of the point k = 0, 
however, for the fermion problem at hand the surface |k| = k-p is relevant, and the cutoff has to 
be defined with respect to this surface. In order to achieve this one proceeds as follows: 

1. Starting from a cutoff-dependent action S'(A) one eliminates all degrees of freedom between 
A and A/s, where s is a factor larger than unity. This gives rise to a new action S"(A' = A/s). 

2. One performs a "scale change" k sk. This brings the cutoff back to its original value and 
a new action S"(A) is obtained. Because of the degrees of freedom integrated out, coupling 
constants (or functions) are changed. 

3. One chooses a value of s infinitesimally close to unity: s = 1 + e, and performs the first two 
steps iteratively. This then gives rise to differential equations for the couplings, which (in 
favorable circumstances) can be integrated until all non-interesting degrees of freedom have 
been eliminated. 



3.1 One dimension 

The one-dimensional case, which has interesting physical applications, will here be mainly used 
to clarify the procedure. Let us first consider a noninteracting problem, e.g. a one-dimensional 
tight-binding model defined by 

^ = Z!^fc4«fc ' ^fc = -2t cos /c - , (3.2) 

k 

where t is the nearest -neighbor hopping integral. We will consider the metallic case, i.e. the 
chemical potential is somewhere in the middle of the band. Concentrating on low-energy prop- 
erties, only states close to the "Fermi points" ±/cf are important, and one can then linearize the 
dispersion relation to obtain 

H= ^F{rk - kp)al^a^^ , (3.3) 

fc,r=it 

where vp = 2tsinkp is the Fermi velocity, and the index r differentiates between right- and left- 
going particles, i.e. particles close to kp and —kp. To simplify subsequent notation, we (i) choose 
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energy units so that vp = 1, (ii) translate /c-space so that zero energy is at A; = 0, and (iii) replace 
the fc-sum by an integral. Then 

/A f]U 
^-rkal{k)aAk) . (3.4) 

For the subsequent renormalization group treatment we have to use a functional integral for- 
mulation of the problem in terms of Grassmann variables (a detailed explanation of this formalism 
is given by Negele and Orland The partition function becomes 



Z(A) = j P0e-^(^) , (3.5) 
where "Dcj) indicates functional integration over a set of Grassmann variables. The action is 

5(A) = f/r|^y_^-0:(A;,r)9.0,(fc,r) + /J(0*,0)| , (3.6) 

where the zero-temperature limit /? — oo has to be taken, and H{(j)*, 0) indicates the Hamiltonian, 
with each replaced by a 0*, and each a replaced by a (p. Fourier transforming with respect to 
the imaginary time variable 

Mk, r)=Tj2 Mk, ^n)e-'""" (cOn = 27r(n + 1/2)T) (3.7) 

and passing to the limit T ^ one obtains the noninteracting action 

1^ / ^€ik,u^)Hu^ + rk]Mk,uj) . (3.8) 
-OO Zn J-A Ztt 

We notice that this is diagonal in k and u which will greatly simplify the subsequent treatment. 
Because of the units chosen, u has units of (length)"^ (which we will abbreviate as L~^), and then 
(f)r{k,uj) has units L^/^. 

We now integrate out degrees of freedom. More precisely, we will integrate over the strip 
A/s < \k\ < A, — OO < uj < OO. The integration over all uj keeps the action local in time. One 
then has 

Z(A) = Z(A,A/s)Z(A/s) , (3.9) 

where Z(A,A/s) contains the contributions from the integrated degrees of freedom, and Z{K/s) 
has the same form of eq. (^.51). The new action is then Sq{K/s) = So{A/s). Introducing the scale 
change 

k' = ks, uj' = ujs, 0' = 0s"3/2 (3 10) 

one easily finds that S'q{A) = 5*0 (A). The action does not change therefore under scale change (or 
renormalization): we are at a fixed point. One should notice that the scale change of k implies 
that k' is quantized in in units of Ak' = 27rs/L, i.e. eliminating degrees of freedom actually 
implies that we are considering a shorter system, with correspondingly less degrees of freedom. 
This means that even though the action is unchanged the new Z{A) is the partition function of a 
shorter system. To derive this in detail, one has to take into account the change in the functional 
integration measure due to the scale change on 0. 
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Before turning to the problem of interactions, it is instructive to consider a quadratic but 
diagonal perturbation of the form 

— -f,{k,u;)<p:{k,u;)<f>,{k,u;) . (3.11) 

We assume that nik^u) can be expanded in a power series 

/i(A;,u;) = /xoo + Atio^ + /ioiit^ + • • • (3.12) 
Under the scale change ( ^.l(j| ) one then has 

f^nm ^ ^ l^nm • (3.13) 

There now are three cases: 

1. a parameter ^nm grows with increasing s. Such a parameter is called relevant. This is the 
case for /loo- 

2. a parameter remains unchanged (/iio,/^oi)- Such a parameter is marginal. 

3. Finally, all other parameter decrease with increasing s. These are called irrelevant. 

Generally, one expects relevant parameters, which grow after elimination of high-energy degrees of 
freedom, to strongly modify the physics of the model. In the present case, the relevant parameter 
is simply a change in chemical potential, which doesn't change the physics much (the same is true 
for the marginal parameters). One can easily see that another relevant perturbation is a term 
coupling right- and left-going particles of the form m(0^02 + 020i)- This term in fact does lead 
to a basic change: it leads to the appearance of a gap in the spectrum. 

Let us now introduce fermion-fermion interactions. The general form of the interaction term 
in the action is 

Si= f m(1234)0*(1)0*(2)0(3)0(4) . (3.14) 

Here 0(3) is an abbreviation for 0^3(^3, "^s), and similarly for the other factors, while u is an 
interaction function to be specified. The integration measure is 

/ = fn /" ^ P] ^(^1 + ^2 - fcs - ^4)5(^1 + 072 - ^3 - W4) . (3.15) 

We now note that the dimension of the integration measure is L~^, and the dimension of the 
product of fields is L^. This in particular means that if we perform a series expansion of u 
in analogy to eq.( 3.12|) the constant term will be s-independent, i.e. marginal, and all other 
terms are irrelevant. In the following we will thus only consider the case of a constant {k- and 
cj-independent) u. 

These considerations are actually only the first step in the analysis: in fact it is quite clear that 
(unlike in the noninteracting case above) integrating out degrees of freedom will not in general 
leave the remaining action invariant. To investigate this effect, we use a more precise form of the 
interaction term: 

Si = £E{^?iC+(1)0:'-(2)0..+ (3)0,_(4)+^720:-,(1)0:,42)0..-(3)0.+ (4)} . (3.16) 
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Here we have reintroduced spin, and the two couphng constants gi and g2 denote, in the original 
language of eg. (^.31), backward ((fcp, — ^f) — ^ {—kpykF)) and forward ((/cp, — ^f) (^f,— ^f)) 
scattering. Note that in the absence of spin the two processes are actually identical. 
Now, the Kadanoff- Wilson type mode elimination can be performed via 

e"^' = J V^e-^ , (3.17) 

where T>(j) denotes integration only over degrees of freedom in the strip A/s < \k\ < A. Dividing 
the field (j) into (p (to be eliminated) and 0' (to be kept), one easily sees that the noninteracting 
action can be written as 5*0 = 5'o(0') + 5*0(0). For the interaction part, things are a bit more 
involved: 

Si = j2Si, = Si,o + Si . (3.18) 

1=0 

Here Sj^i contains i factors (p. We then obtain 

g-s' ^ g-5o{0')-S/,o J v<pe-^°^^^-^' . (3.19) 

Because Si contains up to four factors 0, the integration is not straightforward, and has to be 
done via a perturbative expansion, giving 



j V(j)e-^°^^^~^' = Zo(A,A/s)exp 



OO -| 



*=1^ 



(3.20) 



where the notation (. . .)o,con indicates averaging over and only the connected diagrams are to 
be counted. It can be easily seen, moreover, that because of the U{1) invariance of the original 
action associated to the particle number conservation, terms which involve an odd number of 0' 
or fields are identically zero. The first order cumulants give corrections to the energy and the 
chemical potential and are thus of minor importance. The important contributions come from the 
second order term {S]2)o,con which after averaging leads to terms of the form 0'*0'*0'0', i.e. to 
corrections of the interaction constants gi^2- The calculation is best done diagrammatically, and 
the four intervening diagram are shown in fig.^ 

One can easily see that not all of these diagrams contribute corrections to gi or g2. Specifically, 
one has 

5gi oc gig2[ia) + (c)] + 2gf{d) 

oc {gi+gl){a)+g',{b) (3.21) 



/2 



where the factor 2 for diagram {d) comes from the spin summation over the closed loop. Because 
the only marginal term is the constant in m(1234), one can set all external energies and momenta 
to zero. The integration over the internal lines in diagram (a) then gives 



(a) 



dk r duo 1 1 



2tt J 2tt \uj — k — k 

^ cifcl 1 , , 

— - = —di , 3.22 

A/s 2'Kk 2tt ' ^ ^ 
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where s = 1 + di, and similarly the particle-hole diagrams [b) to [d) give a contribution —di/{27r). 
Performing this procedure recursively, using at each step the renormalized couplings of the previous 
step, one obtains the renormalization group equations 



TT 



(3.23) 



where s = e^. These equations describe the effective coupling constants to be used after degrees 
of freedom between A and Ae^ have been integrated out. As initial conditions one of course 
uses the bare coupling constants appearing in eq. (|3.16|) . Equations (|3.23|) are easily solved. The 



(a) 



(b) 




Figure 4: The diagrams intervening in the renormahzation of the coupHng constants gi and 52 • Note 
that in (b) the direction of one arrow is reversed with respect to (a), i.e. this is a particle-hole diagram. 

combination gi — 2g2 is ^-independent, and one has further 

9ii£) = • (3.24) 

1 + gii 

There then are two cases: 

1. Initially, gi > 0. One then renormalizes to the fixed line gl = gi{i —>■ 00) = 0, g2 = g2 — gi/2, 
i.e. one of the couplings has actually vanished from the problem, but there is still the free 
parameter g2. A case like this, where perturbative corrections lead to irrelevancy, is called 
"marginally irrelevant" . 

2. Initially, gi < 0. Then gi diverges at some finite value of i. We should however notice that, 
well before the divergence, we have left the weak-coupling regime where the perturbative 
calculation leading to the eq. ( p. 23] ) is valid. We should thus not overinterpret the divergence 



and just remember the renormalization towards strong coupling. This type of behavior is 
called "marginally relevant" . 

We will discuss the physics of both cases in the next section. 

Two remarks are in order here: first, had we done a straightforward order-by-order pertur- 
bative calculation, integrals like eq.( p.22| ) would have been logarithmically divergent, both for 
particle-particle and particle-hole diagrams. This would have lead to inextricably complicated 
problem already at the next order. Secondly, for a spinless problem, the factor 2 in the equation 
for gi [£) is replaced by unity. Moreover, in this case only the combination gi — g2 is physically 
meaningful. This combination then remains unrenormalized. 
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3.2 Two and three dimensions 



We will now follow a similar logic as above to consider two and more dimensions. Most argu- 
ments will be made for the two-dimensional case, but the generalization to three dimensions is 
straightforward. The argument is again perturbative, and we thus start with free fermions with 
energy 

^j^ = — - ^ = VFk + 0{k'^) {vF = kF/m). (3.25) 

We use upper case momenta K to denote momenta measured from zero, and lower case to denote 
momenta measured from the Fermi surface: /c = |K| — kp. The Fermi surface geometry now is that 
of a circle as shown in fig.|^. One notices in particular that states are now labeled by two quantum 





Figure 5: Fermi surface geometry in two dimensions. 



numbers which one can take as radial (k) and angular (6). Note that the cutoff is applied around 
the low-energy excitations at |K| — /cp, not around K = 0. The noninteracting action then takes 
the form 



^0 =kp r ^ r f- r ^[ct^\keu){-iu - k)ct^{keu)] 



(3.26) 



One notices that this is just a (continuous) collection of one-dimensional action functional, param- 
eterized by the variable 6. The prefactor kp comes from the two-dimensional integration measure 
(PK = {kp + k)dkd6, where the extra factor k has been neglected because it is irrelevant, as 
discussed in the previous section. 

The general form of the interaction term is the same as in the one-dimensional case 



Si 



Ku 



w(1234)0*(l)0*(2)0(3)0(4) 



(3.27) 



however, the integration measure is quite different because of two-dimensional K-space. Per- 
forming the integration over K4 and a;4 in the two-dimensional analogue of eq. (|3.15|) , the measure 
becomes 



2n) \IjJ-oo 2n Jo 2tt J-a 27r / ^ ' ^ 



I) 



(3.28) 
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Here K4 = Ki + K2 — K3. Now the step function poses a problem because one easily convinces 
oneself that even when Ki^2,3 are on the Fermi surface, in general K4 can be far away from it. 
This is quite different from the one-dimensional case, where everything could be (after a trivial 
transformation) brought back into the vicinity of k = 0. 

To see the implications of this point, it is convenient to replace the sharp cutoff in eq.( |3.28| ) 
by a soft cutoff, imposed by an exponential: 

0(A- |A;4|) ^ exp(-|A:4|/A) . (3.29) 

Introducing now unit vectors fij in the direction of Kj via Kj = {kp + ki)fli one obtains 

ki = \kF{ni + n2- rig) + kiCii + ...{-kp^ fcpdAi - 1) a = rii + ria - rig . (3.30) 

Now, integrating out variables leaves us with A A/s in eq. (|3.28| ) everywhere, including the 
exponential cutoff factor for k^. After the scale change ( p.lOj ) the same form of the action as 
before is recovered, with 

u\ki, 00',, e[) = e-(^-i)('=^/^)ll^l-ilM(fe,/s, ujs, Oi) . (3.31) 

We notice first that nothing has happened to the angular variable, as expected as it parameterizes 
the Fermi surface which is not affected. Secondly, as in the one-dimensional case, the k and uo 
dependence of u is scaled out, i.e. only the values m(0, 0, 6i) on the Fermi surface are of potential 



interest (i.e. marginal). Thirdly, the exponential prefactor in eq. ( |3.31| ) suppresses couplings for 
which |A| 7^ 1. This is the most important difference with the one-dimensional case. 
A first type of solution to |A| = 1 is 

fii = fig ^ = 5 or 

rii = ri4 ^ ria = rig . (3.32) 

These two cases only differ by an exchange of the two outgoing particles, and consequently there 
is a minus sign in the respective matrix element. Both processes depend only on the angle 9i2 
between fii and fia, and we will write 

u(0, 0, ^1, ^2, ^1, ^2) = -m(0, 0, e^, 02, 02, Oi) = FiOi - 62) . (3.33) 

We can now consider the perturbative contributions to the renormalization of F. To lowest 
nontrivial (second) order the relevant diagrams are those of Fermi liquid theory and are reproduced 
in fig.^a. 

Consider diagram (a). To obtain a contribution to the renormalization of F, both Q and 
P — Q have to lie in the annuli to be integrated out. As can be seen from fig.^J, this will give a 
contribution of order di'^ and therefore does not contribute to a renormalization of F. The same 
is true (if we consider the first case in eq.( ]3.32| )) for diagram (b). Finally, for diagram (c), because 



K is small, the poles of both intervening Green functions are on the same side of the real axis, 
and here then the frequency integration gives a zero result. For the second process in eq.( p.32|) the 



same considerations apply, with the roles of diagrams (b) and (c) interchanged. The conclusion 
then is that F is not renormalized and remains marginal: 

f = 0. (3.34) 
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P-Q 



4 



2 




(a) 



(b) 



Figure 6: (a) second order diagrams renormalizing the coupling function. Here P = Ki + K2,K3 = 
Ki — K, and K4 = Ki — K'. Q is the loop integration variable. Note that these diagrams are actually 
identical to those occurring in Landau's theory (figj^). (b) phase space for diagram (a). The rings are 
the degrees of freedom to be integrated out between A and A/s. Note that only if P = are Q and 
P — Q simultaneously in the area to be integrated, giving a contribution of order di. 

The third possibility is to have fli = —0,2, O3 = — ^^4- Then the angle between fii and fl^ 
can be used to parameterize u: 



In this case P = 0, and therefore in diagram (a) if Q is to be eliminated, so is — Q. Conse- 
quently, one has a contribution of order di. For the other two diagrams, one finds again negligible 
contributions of order di'^. Thus, one obtains 



This is a renormalization equation for a function, rather than for a constant, i.e. one here has an 
example of a "functional renormalization group" . Nevertheless, a Fourier transform 




(3.35) 




(3.36) 




(3.37) 



brings this into a more standard form: 



dVx 




di 4tt 
This has the straightforward solution 



Vx 



Vx 



(3.39) 




From eqs.( p.3^ ) and eq.( p.39|) there are now two possibilities: 
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1. At least one of the V\ is negative. Then one has a divergence of Vx{i) at some finite energy 
scale. Given that this equation only receives contributions from BCS-like particle-particle 
diagrams, the interpretation of this as a superconducting pairing instability is straightfor- 
ward. The index A determines the relative angular momentum of the particles involved. 

2. All V\ > 0. Then one has the fixed point Vx = 0, -F(6'i — ^2) arbitrary. What is the underlying 
physics of this fixed point? One notices that here 63 = 61, = 62, i.e. the marginal term in 
the action is (t>*B^4'*e2'^ei^d2- operator language, this translates into 

Hint ~ / deide2ng^ng^ . (3.40) 



We now can recognize this as an operator version of Landau's energy functional, eg. ( |2.13|) . 
The fixed point theory is thus identified as Landau's Fermi liquid theory. 

The generalization of the above to three dimensions is rather straightforward. In addition 
to the forward scattering amplitudes F, scattering where there is an angle 0i2;34 spanned by 
the planes (f2i,f22) and (fl^jfl^) is also marginal. For 0i2;34 7^ these processes are the ones 
contributing to the quasiparticle lifetime, as discussed in sec. |2.2.3| , however they do not affect 
equilibrium properties. The (zero temperature) fixed point properties thus still only depend on 
amplitudes for 0i2;34 = 0, i.e. the Landau /-function. 



4 Bosonization and the Luttinger Liquid 

The Fermi liquid picture described in the preceding two sections is believed to be relevant for 
most three-dimensional itinerant electron systems, ranging from simple metals like sodium to 
heavy-electron materials. The best understood example of non-Fermi liquid properties is that of 
interacting fermions in one dimension. This subject will be discussed in the remainder of these 
lecture notes. We have already started this discussion in section where we used a perturbative 
renormalization group to find the existence of one marginal coupling, the combination gi — 2g2. 
This approach, pioneered by Solyom and collaborators in the early 70's |]l6l, can be extended to 



stronger coupling by going to second or even third order in perturbation theory. A principal 
limitation remains however the reliance on perturbation theory, which excludes the treatment 
of strong-coupling problems. An alternative method, which allows one to handle, to a certain 
extent, strong-interaction problems as well, is provided by the bosonization approach, which will 
be discussed now and which forms the basis of the so-called Luttinger liquid description. It 
should be pointed out, however, that entirely equivalent results can be obtained by many-body 
techniques, at least for the already highly nontrivial case of pure forward scattering . 



4.1 Spinless model: representation of excitations 

The bosonization procedure can be formulated precisely, in the form of operator identities, for 
fermions with a linear energy-momentum relation, as discussed in section To clarify notation, 
we will use a+-(a_-)operators for right-(left-)moving fermions. The linearized noninteracting 
Hamiltonian, eq. (|3.3| ) then becomes 

Ho = VfJ2\ ~ ^F)aV,fca+,fc + {-k - A;F)al,A:'^~,A: > (4-1) 
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and the density of states is N{Ep) = l/{7rvp). In the Luttinger model pO| , pi] , one generahzes 
this kinetic energy by letting the momentum cutoff A tend to infinity. There then are two branches 
of particles, "right movers" and "left movers", both with unconstrained momentum and energy, 
as shown in figure [7[ At least for weak interaction, this addition of extra states far from the Fermi 




Figure 7: Single-particle energy spectrum of the Luttinger model. Occupied states are shown in grey, 
the dark grey area represents the states added to make the model solvable. 

energy is not expected to change the physics much. However, this modification makes the model 
exactly solvable even in the presence of nontrivial and possibly strong interactions. Moreover, 
and most importantly, many of the features of this model carry over even to strongly interacting 
fermions on a lattice. 

We now introduce the Fourier components of the particle density operator for right and left 
movers: 



(4.2) 



The noninteracting Hamiltonian (and a more general model including interactions, see below) can 
be written in terms of these operators in a rather simple form and then be solved exactly. This is 
based on the following facts: 

1. the density fluctuation operators Pa, with a = ±, obey Bose type commutation relations: 

aqL 



[pa{-q),Pa'{q')] = Saa'S, 



27r 



(4.3) 



The relation ( [4.3|) for q ^ q' oi a ^ a' can be derived by straightforward operator algebra. 
The slightly delicate part is eq.( f4.3| ) for q = q' . One easily finds 



[p+(-9),P+(9)] 



k 



nk) 



(4.4) 



where fik is an occupation number operator. In a usual system with a finite interval of states 
between —kp and kp occupied, the summation index of one of the n operators could be 
shifted, giving a zero answer in eq. ([4.4|) . In the present situation, with an infinity of states 
occupied below kp, this is not so. Consider for example the ground state and g > 0. Then 



each term in eq.( [4.4| ) with kp < k < kp + q contributes unity to the sum, all other terms 
vanish, thus establishing the result (O). More generally, consider a state with all levels 
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below a certain value feg (< kp) occupied, but an arbitrary number of particle hole pairs 
excited otherwise. One then has, assuming again q > 0, 

k 




= H nk = ^ . (4.5) 

The result is independent of /cq, and one thus can take the limit /cq ~^ — oo. Together with an 
entirely parallel argument for p_, this then proves eg. ( |4.3| ). Moreover, for g > both p+(— g) 
and p- (g) annihilate the noninteracting groundstate. One can easily recover canonical Bose 



commutation relations by introducing normalized operators, e.g. p+(g) = \j2TT/{qL)p^{q) 
would be a canonical creation operator, but we won't use this type of operators in the 
following. 

2. The noninteracting Hamiltonian obeys a simple commutation relation with the density op- 
erators. For example 

[Ho, pa{q)] = vpaqpaiq) > (4-6) 

i.e. states created by p+(g) are eigenstates of Hq, with energy vpq. Consequently, the kinetic 
part of the Hamiltonian can be re-written as a term bilinear in boson operators, i.e. quartic 
in fermion operators: 

Ho = ^ E P»iQ)Pai-q) . (4.7) 

This equivalence may be made more apparent noting that p+{q) creates particle-hole pairs 
that all have total momentum q. Their energy is Ek+q — £k, which, because of the linearity 
of the spectrum, equals vpq, independently of k. Thus, states created by p+(g) are linear 
combinations of individual electron-hole excitations all with the same energy, and therefore 
are also eigenstates of ([4.1|). 

3. The above point shows that the spectra of the bosonic and fermionic representations of Hq 
are the same. To show complete equivalence, one also has to show that the degeneracies of 
all the levels are identical. This can be achieved calculating the partition function in the 
two representations and demonstrating that they are equal. This then shows that the states 
created by repeated application of p± on the ground state form a complete set of basis states 



We now introduce interactions between the fermions. As long as only forward scattering of 
the type (fcp; — ^f) — (^f; — ^f) or (fcp; ^f) ~^ (^f; ^f) is introduced, the model remains exactly 
solvable. The interaction Hamiltonian describing these processes takes the form 

i^int = ^ E i92{q)Pa{q)p-a{-q)+g4{q)Pa{q)Pa{-q)j ■ (4.8) 
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Here, 5^2(9) and (74(5) are the Fourier transforms of a real space interaction potential, and in a 
realistic case one would of course have g2{q) = gi{<i) = dil), but it is useful to allow for differences 
between (72 and g^. For Coulomb interactions one expects (?2,5'4 > 0. In principle, the long-range 
part of the Coulomb repulsion leads to a singular g-dependence. Such singularities in the gi can be 
handled rather straightforwardly and can lead to interesting physical effects as will be discussed 
below. Here I shall limit myself to nonsingular 5'2,5'4- Electro n-phonon interactions can lead to 
effectively attractive interactions between electrons, and therefore in the following I will not make 
any restrictive assumptions about the sign of the constants. One should however notice that a 
proper treatment of the phonon dynamics and of the resulting retardation effects requires more 
care 
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Putting together ( |4.7] ) and ( |4.8| ), the complete interacting Hamiltonian, the Tomonaga-Luttinger 
model, then becomes a bilinear form in boson operators that is easily diagonalized by a Bogolyubov 
transformation. A first consequence is the expression for the excitation spectrum 



uj{q) = \q\ 



The diagonal boson operators are linear combinations of the original p operators, and consequently, 
these elementary excitations are collective density oscillations, their energy being determined both 
by the kinetic energy term and the interactions. 

We note here that in order for the Bogolyubov transformation to be a well-defined unitary 
transformation, g2{q) has to decrease at large q at least as |g|~^/^.||23| On the other hand, the large- 



q behavior of g2 is unimportant for the low-energy properties of the model. We therefore in the 
following will almost always use a g-independent g2 and g^. An approximate and frequently used 
way to cure the divergences arising due to this procedure is to keep a parameter a in subsequent 
formulae as a finite short-distance cutoff, of the order of a lattice spacing. One can then also 
include the "backward scattering" (kp; —kp) {—kp', kp), because for spinless electron this is just 
the exchange analogue of forward scattering and does not constitute a new type of interaction. 
It is worthwhile emphasizing here that the solution is valid for arbitrarily strong interactions, no 
perturbative expansion is needed! 

Up to this point, the construction does not allow for a direct calculation of correlation functions 
like the one-particle Green function or more generally any function involving individual creation or 
destruction operators. This type of correlation function becomes tractable by representing single 
particle operators in terms of the boson operators. To this end, we introduce the field operators 



H(x) = ^Ee-"''''/'-^''^(^P+(?)-P-(?)j +^ (4.11) 

N = N+ + N_ , J = N+-N_ . (4.12) 

Here N± are the numbers of particles added to the ground state on the right- and left-moving 
branch i.e. 

N^ = p^{q = 0)=J2 in±,k - {n±,k)o) • (4.13) 

k 
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Because addition of a particle changes both and J, one has the "selection rule" (—1)^ = (—1)"^; 
a is a cutoff parameter which (at least in principle, see the discussion above) has to be set to zero in 
the end of any calculation. The fields and 11 then obey canonical boson commutation relations: 

[0(x),n(y)] = i5(x-y) , (4.14) 

and (f) is related to the local particle density via 

= -vr(p(a;) - po) , (4.15) 

where po is the average particle density in the ground state. More precisely, in a lattice model 
this would represent the slowly varying components (g ~ 0) of the density, whereas components 
with q ^ 2kp would correspond to crossproducts between 'ip±. 
The expression for the single fermion operators then is 

^±(a;) = lim U± exp \±ikpx =F i0(x) + i9{x)] , (4.16) 

where the upper and lower sign refer to right- and left-moving electrons respectively, and 

6(x) = Tc J n(x') dx' 

The ^/-operators (sometimes referred to as "Klein factors"; they are non-Hermitian: ^ Ua) 
decrease the total particle number on one of the branches by unity and are necessary because the 
boson fields all conserve the total particle number. These operators also insure proper anticom- 
mutation between right- and left-going operators: one easily checks that the "chiral components" 
0+ = 6* — and 0_ = + commute with each other, one therefore needs the anticommutation 
relations [Ua, U(^]+ = 0, etc. In the thermodynamic limit L oo, the fact that the f/'s change 
the particle number is of minor importance because this represents a shift of fci? by a quantity of 
order 1/L, and one then can replace the U^s by Majorana (Hermitian) fermion operators obeying 



[Va,Vi3]+ = &s discussed in more detail in the appendix. In single-chain problems these 

effects play a minor role, however in the many chain systems to be discussed at the end of these 
notes, proper account of anticommutation is crucial. 



A detailed derivation of the important eq. ( [4.16| ) as an operator identity is given in the literature 



23| , |2^ . However, a simple plausibility argument can be given: creating a particle at site x 
requires introducing a kink of height vr in 0, i.e. has to be shifted by vr at points on the left of 
X. Displacement operators are exponentials of momentum operators, and therefore a first guess 
would be '0(a;) ~ exp(i7r /f^ n(x')(ix'). However, this operator commutes with itself, instead of 
satisfying canonical anticommutation relations. Anticommutation is achieved by multiplying with 
an operator, acting at site x, that changes sign each time a particle passes through x. Such an 



operator is exp(±i0(x)). The product of these two factors then produces (|4.16|) . 

The full Hamiltonian can also be simply expressed in terms of and H. In the long-wavelength 
limit, neglecting the momentum dependence of the Qi, one can express the total Hamiltonian in 
the field phase formalism 

H = Ho + ifint = I dx\ ^U{xy + ^id.<P? 1 ■ (4.18) 
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This is obviously just the Hamiltonian of an elastic string, with the eigenmodes corresponding 
to the collective density fluctuations of the fermion liquid. It is important to notice that these 
collective modes are the only (low-energy) excited states, and that in particular there are no 
well-defined single particle excitations, nor are there the incoherent particle-hole pair excitations 
typical of a Fermi gas. The parameters in ( ^.18| ) are given by 



u = JU + - (^Y . K = J T"'^''-'' . (4,19) 



2tt J \2tx J ^ ^txvy + gi + §2 

The energies of the eigenstates are uj{q) = u\q\, in agreement with eq. ( [4.9| ). From the continuity 
equation, the expression (|4.15|) for the local particle density and the equation of motion of the 
(number) current density is 

]{x) = uKU{x) . (4.20) 

Note in particular that for g2 = g^ one has uK = vp, i.e. the expression for the current density is 
interaction-independent. The relation uK = vp holds in particular for systems with full (Galilean) 
translational invariance. On the other hand, in the continuum limit of lattice systems this relation 
is in general not true. 

The most remarkable result here is the "collectivization" of the dynamics: there are no 
quasiparticle-like excitations. In fact there is a rather simple physical picture explaining this: 
imagine accelerating one particle a little bit in one direction. Very soon it will hit its neighbor 
and transmit its momentum to it, and the neighbor will in turn transmit its momentum to a 
further neighbor, and so on. Quite quickly, the initial localized motion will have spread coherently 
through the whole system. This picture can be formalized noting that in one dimension the differ- 
ence between a solid and a fluid is not well-defined: whereas is higher dimensions solids and fluids 
are differentiated by the presence or absence of long-wavelength transverse modes, no transverse 
modes can exist in a system with movement along only one direction. The long-wavelength modes 
thus can equally well be considered as the phonons of a one-dimensional crystal |2^, Note 
that on the contrary in dimensions larger than one the neighbors of any given particle can be 
pushed aside, giving rise to a backflow that allows the particle to move trough the system more 
or less freely. 

The exclusive existence of collective excitations, at least at low energies, is one of the typical 
properties of the Luttinger liquid. Rather than discussing the physics of the spinless case in detail, 
we will turn now to the more interesting case of fermions with spin. 

4.2 Model with spin; the concept of the Luttinger Liquid 

In the case of spin-1/2 fermions, all the fermion operators acquire an additional spin index s. 
Following the same logic as above, the kinetic energy then takes the form 



^0 = v^Y^ii^- ^F)a+,fc,s«+,fc,s + {-k - A;F)al,fc,,a_,fc,s 

k^s y J 



= ^ E Pa,.(g)p«,.(-g) , (4.21) 

where density operators for spin projections s =t, [ have been introduced: 

P±,^(9) = E4,fc+g,««±,M- (4-22) 

k 
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There are now two types of interaction. First, the "backward scattering" (/cp, s; —kp, t) —>■ 
(— fcp, s; kp, t) which for s ^ t cannot be re-written as an effective forward scattering (contrary to 
the spinless case). The corresponding Hamihonian is 

-f^int,! = X! 9l(^+,k,s^-,p,t(^+,p+2kp+q,t^-,k-2kF-q,s ■ (4.23) 

k,p,q,s,t 

And, of course, there is also the forward scattering, of a form similar to the spinless case 

^mt,2 = 177^2 ( ^2(g)Pa,s(g)p-a,t(-g) +^4(g)Pa,s(g)Pa,t(-g) ) ■ (4.24) 

^-^ q,a,s,t \ / 

To go to the bosonic description, one introduces and 11 fields for the two spin projections 
separately, and then transforms to charge and spin bosons via (pp^a- = (0t 4>l)/V^^ ^p,a = 
(Il-f ± nj/v^. The operators (p^ and III, obey Bose-like commutation relations: 

[Mx),U^{y)] = i6,,6{x-y) , (4.25) 

and single fermion operators can be written in a form analogous to ( [4.161 ): 



ij±,s{x) = lim exp I iifcpx - i{±{(j)p + s(j)a) + {Op + sO^) )/V2 1 , (4.26) 



-0 
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where 9^{x) = n Il^{x')dx' . 

The full Hamiltonian H = Hq + Hint^i + Hint,2 then takes the form 

H = Hp + H, + J dx cos( V80,) . (4.27) 

Here a is a short-distance cutoff, and for u = p,(T 

H.= ldx\ ^-^lil + -^{dA.f 1 , (4.28) 



with 



and Qp = gi — 2(72, Qa = Qi, g4,p = 94, gi,a = 0. The choice of sign (which is the conventional one) 
for the cosine-term in eq. (|4.27|) corresponds to a particular "gauge choice", as discussed in the 



appendix. 

For a noninteracting system one thus has Uy = v-p (charge and spin velocities equal!) and 
Ky = 1. For Qi = 0, ([4.27|) describes independent long- wavelength oscillations of the charge and 
spin density, with linear dispersion relation uj^ik) = Uu\k\, and the system is conducting. As 
in the spinless case, there are no single-particle or single particle-hole pair excited states. This 
model (no backscattering) , usually called the Tomonaga-Luttinger model, is the one to which the 
bosonization method was originally applied [^, |29[| . 

For gi ^ the cosine term has to be treated perturbatively. We have already obtained the 
corresponding renormalization group equations in the previous section (eq.( |3.23D ). In particular, 
for repulsive interactions {gi > 0), gi is renormalized to zero in the long- wavelength limit, and 
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at the fixed point one has K^. = 1. The three remaining parameters in ([4.27| ) then completely 
determine the long-distance and low-energy properties of the system. 

It should be emphasized that ([4.27|) has been derived here for fermions with linear energy- 
momentum relation. For more general (e.g. lattice) models, there are additional operators 
arising from band curvature and the absence of high-energy single-particle states [^. One can 
however show that all these effects are, at least for not very strong interaction, irrelevant in the 
renormalization group sense, i.e. they do not affect the low-energy physics. Thus, l \i-2T{ ) is still 



the correct effective Hamiltonian for low-energy excitations. The lattice effects however intervene 
to give rise to "higher harmonics" in the expression for the single-fermion operators, i.e. there are 



low energy contributions at wavenumbers q ~ (2m + l)kY for arbitrary integer m ||27 



The Hamiltonian ( [4.27| ) also provides an explanation for the physics of the case of negative 



(y^i, where the renormalization group scales to strong coupling (eq. (|3.23| )). In fact, if \gi\ is large 
in (|4.27|) , it is quite clear that the elementary excitations of 00- will be small oscillations around 



one of the minima of the cosine term, or possibly soliton-like objects where 0o- goes from one of 
the minima to the other. Both types of excitations have a gap, i.e. for < one has a gap in 
the spin excitation spectrum, whereas the charge excitations remain massless. This can actually 



investigated in more detail in an exactly solvable case |Q. We will subsequently concentrate 
on the case gi > 0, so that there is no spin gap, however investigations of spectral functions as 
described below have also been recently performed for the case with a spin gap. [pl[] 

4.2.1 Spin— charge separation 

One of the more spectacular consequences of the Hamiltonian ([4.27|) is the complete separation of 
the dynamics of the spin and charge degrees of freedom. For example, in general one has ^ Up, 
i.e. the charge and spin oscillations propagate with different velocities. Only in a noninteracting 
system or if some accidental degeneracy occurs one does have u„ = Up = vp. To make the meaning 
of this fact more transparent, let us create an extra particle in the ground state, at t = and spatial 
coordinate Xq. The charge and spin densities then are easily found, using p(x) = —{\^/n)dx(f)p 
(note that p{x) is the deviation of the density from its average value) and az{x) = —{\^/Tr)dx(j)a 

{0\ilj+{xo)pix)ip'l{xo)\0) = 6{x-Xo) , 
(0|V^+(xo)a,(x)^i(xo)|0) = 5(x-xo) . (4.30) 

Now, consider the time development of the charge and spin distributions. The time-dependence 
of the charge and spin density operators is easily obtained from (|4.271 ) (using the fixed point value 
gi =0), and one obtains 

{0\tp+{xo)p{x,t)iji{xo)\0) = 5{x-XQ-Upt) , 

{Q\^+{xo)a,{x,t)ijl{x^)\Q) = 6{x-xo-uJ) . (4.31) 

Because in general Mo- 7^ Up, after some time charge and spin will be localized at completely 
different points in space, i.e. charge and spin have separated completely. 

Here a linear energy-momentum relation has been assumed for the electrons, and consequently 
the shape of the charge and spin distributions is time-independent. If the energy-momentum 
relation has some curvature (as is necessarily the case in lattice systems) the distributions will 
widen with time. However this widening is proportional to and therefore much smaller than 
the distance between charge and spin. Thus, the qualitative picture of spin-charge separation is 
unchanged. 
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4.2.2 Physical properties 



The simple form of the Hamiltonian ( [4.27| ) at the fixed point gl = makes the calculation of 
physical properties rather straightforward. The specific heat now is determined both by the 
charge and spin modes, and consequently the specific heat coefficient 7 is given by 



7/70 = 7t(^'f/"p + ^p/^ia) 



(4.32) 



Here 70 is the specific heat coefficient of noninteracting electrons of Fermi velocity vp. 

The spin susceptibility and the compressibility are equally easy to obtain. Note that in ( 4.27| ) 
the coefficient u^j/K^ determines the energy necessary to create a nonzero spin polarization, and, 
as in the spinless case, Up/Kp fixes the energy needed to change the particle density. Given the 
fixed point value K* = 1, one finds 



X/XO = Vp/Ua , l^/l^Q = VpKp/Up 



(4.33) 



where Xo '^o are the susceptibility and compressibility of the noninteracting case. From 
eqs.( |4.32|) and ([4.33|) the Wilson ratio is 



p _ A /U 

-rtw — 



XTo 
1X0 



2Ur 



Up + u„ 



(4.34) 



The quantity Tlp{x) is proportional to the current density. As before, the Hamiltonian com- 



mutes with the total current, one thus has 



criuj) = 2KpUp5{ijj) + a 



rcg 



(4.35) 



i.e. the product KpUp determines the weight of the DC peak in the conductivity. If the total 
current commutes with the Hamiltonian cxreg vanishes, however more generally this part of the 
conductivity varies as cu^ at low frequencies |^ . 



The above properties, linear specific heat, finite spin susceptibility, and DC conductivity are 
those of an ordinary Fermi liquid, the coefficients Up,u^, and Kp determining renormalizations 
with respect to noninteracting quantities. However, the present system is not a Fermi liquid. This 
is in fact already obvious from the preceding discussion on charge-spin separation, and can be 
made more precise considering the single-particle Green function. Using the representation (|4.26| ) 
of fermion operators one finds (at the fixed point gi = 0) 



G''{x,t) = -ie(t)( V'+,.(a;,t),<s(0,0) 



Re 



1 



TT 



{upt — x){u^t — x) 



a 



[a 



lUptf 



x^ 



5/2\ 



(4.36) 



where the presence of the Heaviside function ensures the retarded nature of and [, ]+ denotes 
the anticommutation of fermion operators. Note that this expression factorizes into a spin and a 
charge contribution which propagate with different velocities. Fourier transforming ( ^.36| ) gives 
the momentum distribution function in the vicinity of k^: 



nk ~ — const, x sign(A; — kF)\k — kp] 



(4.37) 
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and for tlie single-particle density of states (i.e. the momentum-integrated spectral density) one 
finds: 



N{uj) 



\UJ\ 



(4.38) 



In both cases 6 = {Kp + 1/Kp — 2)/4. Note that for any Kp ^ 1, i.e. for any nonvanishing interac- 
tion, the momentum distribution function and the density of states have power-law singularities 
at the Fermi level, with a vanishing single particle density of states at E-p. This behavior is obvi- 
ously quite different from a standard Fermi liquid which would have a finite density of states and 
a step-like singularity in n^. The absence of a step at kp in the momentum distribution function 
implies the absence of a quasiparticle pole in the one-particle Green function. In fact, a direct 
calculation of the spectral function A{k,Lj) from (|4.36|) ||33| , |3^ shows that the usual quasiparti- 
cle pole is replaced by a continuum, with a lower threshold at mm,y{uu{k — kp)) and branch cut 
singularities at = Upp and uj = u^p: 



{uj - u^{k - kp)) 



5-1/2 



UJ 



uJk 



A{k,uj) 

A{k,u) ^{u-Up{k-kp))^^-^^/^ LJ~U^{k-kp) 



kp] 



(5-l)/2 
S-1/2 



{Up > Ua) 
{Up < U„) . 



(4.39) 
(4.40) 



The coefficient Kp also determines the long-distance decay of all other correlation functions of 
the system: Using the representation (|4.26|) the charge and spin density operators at 2kp are 



OcDV\[{x) 



5^^1,s(x)V^. 

s 

lim 



(x) 



i%/2,j>p{x) 



s 

lim e~' 



cos[a/20<^(x)] 



^,-s{x) 

^■^"("^ cos[v^e,(a;)] 



(4.41) 



(4.42) 



Similar relations are also found for other operators. It is important to note here that all these 
operators decompose into a product of one factor depending on the charge variable only by another 
factor depending only on the spin field. Using the Hamiltonian ( [4.27| ) at the fixed point gl = Q 
one finds for example for the charge and spin correlation functions^ 



(n(x)n(O)) = Kp/{-nxf + Aicos{2kpx)x'^~'^nn-^''^{x) 
+ A2 cos{Akpx)x~'^^'' + ... , 
(S(a;)-S(0)) = l/^-Kx)^ + Bicofi{2kpx)x'^-^nYi^'^{x) + ... 



(4.43) 
(4.44) 



with model dependent constants Ai,Bi. The ellipses in ( 4.43 ) and ( 4.44| ) indicate higher harmonics 
of cos(2/cfx) which are present but decay faster than the terms shown here. Similarly, correlation 
functions for singlet (SS) and triplet (TS) superconducting pairing are 



(OL(x)Oss(O)) 
(O^sJx)Ots.(O)) 



Cx-i-i/^Mn-^/^i 



x] 



x) + . 



(4.45) 



^The time- and temperature dependence is also easily obtained, see ref.[p5[ 
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The logarithmic corrections in these functions |^ have been studied in detail recently ^ 



P| . The corresponding susceptibilities (i.e. the Fourier transforms of the above correlation 
functions) behave at low temperatures as 

XcDw(T) ^ T^"-^! HT)r'/' , xsDw(T) ^ T^^'^l HT)\'/' , (4.46) 
Xss(T)^Ti/^''-i|ln(r)|-3/2 , XTs{T)^T'/^^''\HT)\'/\ (4.47) 

i.e. for Kp < 1 (spin or charge) density fluctuations at 2kp are enhanced and diverge at low 
temperatures, whereas for Kp > 1 pairing fluctuations dominate. The "phase diagram" , showing 
in which part of parameter space which type of correlation diverges for T — > is shown in fig.§ 





gl 






SDW 




TS 




(CDW) 


(SS) 




1/2 


2 K 


CDW 


CDW 


SS 


SS 




(SS) 


(CDW) 





Figure 8: Phase diagram for interacting spin~l/2 fermions. 

These correlation functions with their power law variations actually determine experimentally 
accessible quantities: the 2kp and 4kp charge correlations lead to X-ray scattering intensities 
hkj, ~ T^p, I^kp ~ T^-'^p"^, and similarly the NMR relaxation rate due to 2kp spin fluctuations 
varies as 1/Ti fa T^^. The remarkable fact in all the above results is that there is only one 
coefficient, Kp, which determines all the asymptotic power laws. 

We here re-emphasize the two important properties of spin-1/2 interacting fermions in one 
dimension: (i) correlation functions show power-law decay, with interaction-dependent powers 
determined by one coefficient, Kp] and (ii) "spin-charge separation": spin and charge degrees of 
freedom propagate with different velocities. Both these properties are typical of the Luttinger 
liquid and invalidate the Landau quasiparticle concept in one dimension. 

A nice experimental example of different spin-charge separation (different velocities for the 
spin and charge modes) is provided by Raman scattering experiments on single-channel quantum 



wires, pi] , |4^ , ^ On the other hand, the situation in quasi-one-dimensional conductors [0] is less 



clear: in compounds of the TTF-TCNQ series, the observation of strong Akp CDW fluctuations 



45| seems only to be explainable in the Luttinger-liquid picture. On the other hand, in the much 
studies (TMTSF)2X family of compounds, the NMR data can be interpreted in the Luttinger 
liquid framework |H6[, but much of the magnetotransport work of Chaikin and collaborators is 



explained rather satisfactorily using a Fermi liquid like picture. 
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4.2.3 Long— range interactions: Wigner crystallization 

The above calculations can be straightforwardly generalized to the case of long-range interactions. 
Of interest is the case of unscreened Coulomb interactions {V{r) = e^/r) which for example is of 
relevance for the physics of an isolated quantum wire. The short-distance singularity has to be cut 
off, and for example in a wire of diameter d an approximate form would be V{r) = e"^ / \/x^~+~d^ , 
leading to a Fourier transform V{q) = 2e'^Ko{qd). The long-range nature of the interaction is 
only of importance for the forward-scattering processes, and these appear only in the charge part 
of the Hamiltonian which is now given by 



2 



+ \ f dxdx'V{x - x')dr,(j)pd^'(f)p , (4.48) 

where gi = gi/{27[Vp). The elementary excitations then are found to be charge oscillations (plas- 
mons), with energy-momentum relation 



^p(g) = ^F|g| V (1 + 9i){l -~9i + 2V{q) ) (4.49) 

where V{q) = V{q)/{7('Vp). The long-wavelength form (g 0), iOp{q) ~ In gl"*^/^, agrees with 
RPA calculations however, the effect of gi, which is a short-range exchange contribution, 

are usually neglected in those calculations. The spin modes are still given by ujf^{q) = Ua-\q\, with 



Ua = vy\Ii - gl- 

In the evaluation of correlation functions, the charge averages lead to 



m^)-um = I dq^ '^'^"^ 







9i 



q ^ 1 - ^1 + 2V{q) 
!=s C2Vlnx , (4.50) 

with C2 = ^(1 + gi)7ivp/e'^. One thus obtains for example 

(p(x)/o(0)) = Ai cos{2kpx) exp(— C2Vln x)/x 

+A2Cos(4A;Fa;)exp(-4c2yina;) + ... , (4-51) 

where 2 are interaction dependent constants, and only the most slowly decaying Fourier com- 
ponents are exhibited. The most interesting point here is the extremely slow decay (much slower 
than any power law!) of the 4kp component, showing an incipient charge density wave at wavevec- 
tor Akp. This slow decay should be compared with the power-law decay found for short-range 
interactions (eg. ( [4.43] )). The 4kp oscillation period is exactly the average interparticle spacing, 
i.e. the structure is that expected for a one-dimensional Wigner crystal. Of course, because of 
the one-dimensional nature of the model, there is no true long-range order, however, the ex- 
tremely slow decay of the 4kp oscillation would produce strong quasi-Bragg peaks in a scattering 
experiment. It is worthwhile to point out that this Akp contribution arises even if the Coulomb 
interaction is extremely weak and depends only on the long-range character of the interaction. 
On the other hand, any 2kp scattering is considerably weaker, due to the 1/x prefactor in ( |4.51| ) 
which has its origin in the contribution of spin fluctuations. The 2kp spin correlations equally 
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contain a 1/x factor and thus have the same asymptotic decay as the 2kp charge correlations, 
eg. ( |4.51D . On the other hand, correlation functions that involve operators changing the total 



number of particles (e.g. the single particle Green function and pairing correlation functions) de- 
cay like exp[— est. x (Inx)"^/^], i.e. /asier than any power law. This in particular means that the 
momentum distribution function and all its derivatives are continuous at kp, and there is only 
an essential singularity at kp. 

It is instructive to compare the above result ( |4.51| ), obtained in the limit of weak Coulomb 
interactions, with the case of strong repulsion (or, equivalently, heavy particles). The configuration 
of minimum potential energy is one of a chain of equidistant particles with lattice constant a, and 
quantum effects are expected to lead only to small oscillations in the distances between particles. 
The Hamiltonian then is 

^ = E 1^ + ^ E V"ima)iu, - ui^^r , (4.52) 

I l^m 

where ui is the deviation of particle / from its equilibrium position. In the long- wavelength limit, 
the oscillation of this lattice have energy uj{q) = J2/{ma)eq\ ln(ga)|^/^. The most slowly decaying 



part of the density-density correlation function then is 



(p(x)p(O)) ^ cos(27ra:/a) exp (-2^^2^)1/2 • (4-53) 



Noticing that kp = n / (2a), one observes that the results ( |4.51| ) and ( |4.53| ) are (for gi = 0) identical 



as far as the long-distance asymptotics are concerned, including the constants in the exponentials. 
Eq. ( [4. 511 ) was obtained in the weak interaction limit, whereas ( ^.53) ) applies for strong Coulomb 
forces. Similarly, the small-g limit of the charge excitation energies is identical. 

4.3 Summary 

In this section we have developed the basic bosonization formalism for one-dimensional inter- 
acting fermions and seen some elementary and direct applications to the calculation of some 
physical properties, in particular correlation functions. We have seen that the properties of the 
one-dimensional interacting system, the Luttinger liquid, are fundamentally different from two- 
or three-dimensional Fermi liquids. In particular the elementary excitations are not quasiparti- 
cles but rather collective oscillations of the charge and spin density, propagating coherently, but 
in general at different velocities. This gives rise to the interesting phenomenon of spin-charge 
separation. Finally, and again contrary to the Fermi-liquid case, most correlation functions show 
non-universal powerlaws, with interaction-dependent exponents. However, all these exponents 
depend only on one parameter, Kp, the spin analogue of which, K^^, being fixed to unity by 
spin-rotation invariance [K^^ 7^ 1 is possible if spin rotation invariance is broken). Beyond Kp, 
the only parameters that intervene in the low-energy physics of a Luttinger liquid are the veloc- 
ities of the spin and charge modes, Up^^. In the spinless case only two parameters, K and u are 
involved. Finally, we have seen that long-range (Coulomb) interactions can profoundly modify 
these properties. 

For lack of time and space, we have not touched here upon a number of interesting uses and 
generalizations of the bosonization method, in particular the Kondo effect [4^ and applications 



to Fermi systems in more than one dimension. ^ An extensive recent review on Luttinger 



liquids and bosonization has been given by Voit . ||52|| 
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5 Applications 



We will subsequently discuss some results from applications of the formalism developed to interest- 
ing physical questions: transport and the effect of disorder, and the physics of antiferromagnetic 
spin chains. Here, we will come across the two cases where Luttinger liquid behavior is probably 
best established experimentally: the physics of quantum Hall edge states, upon which we will 
only touch briefly, and quantum spin systems, in particular the spin-1/2 antiferromagnet, which 
we will discuss in some detail, with some reference to experiment. 

5.1 Transport 

5.1.1 Conductivity and conductance 

In the previous section we were concerned with equilibrium properties and correlation functions, 
in order to characterize the different phases possible in a one-dimensional system of interacting 
fermions. Here, we will investigate transport, in particular the DC conductivity. Finite-frequency 
effects have also been investigated, and the reader is referred to the literature [^, ^ . 



To clarify some of the basic notions, let us first consider a Luttinger model in the presence of 
a weak space- and time-dependent external potential (p. The interaction of the fermions with ip 
is described by the extra term 

-f^ext = -e J dxp{x)(p{x,t) (5.1) 

in the total Hamiltonian. We will assume that the external field is slowly varying in space, so that 
in the particle-density operator p only products of either two right- or two left-going fermion 
operators appear but no cross terms. Standard linear response theory tells us that the current 
induced by the potential is given by 



^2 



j[x, t) = f dt' ( dx'Djp{x -x',t- t')^{x', t') , (5.2) 

tl J — CO J 

where the retarded current-density correlation function is given by 
D^,ix,t) = -ie(t)([j(x,t),p(0,0)]) 

- ""'^'-Qit) { 6'ix-u,t) + 5'ix + Upt)\ . (5.3) 
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The second line is the result for spin-1/2 electrons. For spinless fermions one has to make the 
replacement UpKp — > uK/2. 

Let us now first consider the situation where we adiabatically switch on a potential of frequency 
u and wavenumber q along the whole length of the system. From eq.( |5.3| ) one then straightfor- 
wardly obtains the g- and tu-dependent conductivity as 

. N i(c^ + iO+) ,^ ^, 

<y[q,^) = ^UpKp- — — — —. 5.4 

In particular, the real part of the conductivity for constant applied field is 

HeaiXijUj) = ——UpKpdiuj) , (5.5) 
n 
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in agreement with eg. ( |4.35| ) (where units with e'^ = h = 1 were used). 



Applying on the other hand a static field over a finite part of the sample, one obtains a current 
j = 2e^KpU/h, where U is the applied tension. The conductance thus is 

G = —K, , (5.6) 

and depends on Kp only, not on Up. For the noninteracting case Kp = 1 this is Landauer's well- 
known result [Q. Note that interactions affect the value of the conductance. The conductance 
here is independent on the length over which the field is applied. Noting that in dimension d the 
conductance is related to the DC conductivity via G = L'^^'^a, a length-independent conductance 
implies an infinite conductivity in one dimension, in agreement with eq.( |5.3| ). The fact that Up 
does not appear in the expression for G can be understood noting that applying a static field 
over a finite (but large) part of the sample, one is essentially studying the wavenumber-dependent 
conductivity at strictly zero frequency, which from eq.( |5.20|) is given by cr(g, 0) = 2e^Kp6{q)/h, 
indeed independent on Up. On the other hand, applying a field of finite frequency over a finite 
length i, one can see that one measures the conductivity a{q ^ 0, cu) if > Up/i. 

The situation of a finite static potential drop over only a finite part of a wire is clearly difficult 
to realize experimentally. However, the result (p.6|), or more precisely its spinless analogue G = 



Ke'^/h, applies to the chiral Luttinger liquid as realized on the edge of quantum Hall systems. ||55 



p^ , p7[ On the other hand, for a quantum wire connected to measuring contacts which impose a 



potential difference over the whole length of the wire, one obtains ||58| , |59| , |56| , pO 



26^ 

G=-. (5.7) 

independent of Kp, i.e. the momentum-conserving interactions play no role in the dc conductance, 
as to be expected. On the other hand, the value of Kp plays an important role for the effects 



of random and contact potentials. lO, 32| These results permit a consistent explanation of some 



experiments on quantum wires [^, but leave open questions in other cases. [l64| 



5.1.2 Persistent current 

The Luttinger model description can be used straightforwardly to obtain the current induced in 
a strictly one-dimensional ring threaded by a magnetic flux $ [^]. The argument can in fact 



be made very simply: in the one-dimensional geometry, the vector fleld can be removed entirely 
from the Hamiltonian via a gauge transformation, which then leads to the boundary condition 
i/^{L) = exp(27ri$/$o)^(0) for the fermion fleld operator. Here L is the perimeter of the ring, and 
$0 = hc/e. For spinless fermions, this is achieved by replacing 

2$ 

U(x) ^ n(x) + -— (5.8) 



in the bosonization formula, eq.( [4.16D . The total J-dependent part of the Hamiltonian then 
becomes 

TT?/ K 

Hj = ^{J + 2^/%f , (5.9) 



giving rise to a number current 
%dE _uK / 2' 
2^M ~~\ ^ $0 



$0^^ uK f 2$\ 
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At equilibrium, J is chosen so as to minimize the energy. Given that at constant total particle 
number J can only change by two units, one easily sees that the equilibrium (persistent) current 
has periodicity $0, and reaches is maximum value uK/L at $ = $o/2, giving rise to the familiar 
sawtooth curve for the current function of flux. 

For fermions with spin, as long as there is no spin gap {gi > 0), the above results can be taken 
over, with the simple replacement uK 2upKp, the factor 2 coming from the spin degeneracy. 
Note in particular that the persistent current, an equilibrium property, is given by the same 
combination of parameters as the Drude weight in the conductivity. This is an illustration of 
Kohn's result relating the Drude weight to the effect of a magnetic flux through a ring. 

In the case of negative gi, electrons can be transfered from the right to the left-going branch 
only by pairs. Consequently, the periodicity of the current and the ground state energy is doubled 
to 2$05 and the maximum current is equally doubled. This behavior has actually been found in 



numerical calculations 67, 



5.1.3 Quantum Hall edge states 

The Luttinger liquid picture has an interesting application to the physics of the fractional quantum 
Hall effect, as discovered and discussed by Wen |70|. To see how this comes about, consider 
the states available in the different Landau levels in the vicinity of the edge of the quantum Hall 
device, as shown in figure ^ It is clear that low-energy states only exist at the edge (the 




r <x angular momentum 



Figure 9: Quantum states in the lowest Landau level in the vicinity of the edge of a quantum hall device. 
The spatial variation of the confining potential is assumed to be slow on the scale of the magnetic length, 
so that the energies of the different quantum states are determined by the local value of the confining 
potential. For a circular device the angular momentum of a state increases proportionally to its distance 
r from the origin. 

bulk quantum Hall state is well-known to be characterized by a finite excitation gap), and close 
to the Fermi energy (i.e. the edge) the states have a linear dispersion relation. This can be made 
particularly clear if one assumes a disk-shaped sample: the states have a wavefunction 

^z'^^ie'^)'' , (5.11) 

with k increasing linearly with radial position. The angular momentum quantum number k thus 
plays a role very similar to linear momentum in the linear geometry we have assumed up to now. 
One thus can linearize the the dispersion in figure |^ and obtains essentially the spinless model 
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discussed in section |4.1| , the only difference being tliat liere only right-going particles exist. This 
difference is the origin of the term chiral Luttinger liquid (in fact, the left going branch is to be 
found on the opposite edge of the device) . Because there are no left-going particles (or at least they 
can be thought of as being at a macroscopic distance), there also is no right-left interaction, and 
consequently one expects the noninteracting value K = 1. Moreover all the left-going components 
of the fields have to be projected out, for example one has to replace — 0+ = (0 — 6')/2. 

However, straightforward adoption of this scheme leads to trouble: from the preceding subsec- 
tion we know that K = 1 leads to a conductance (which in the present case is the Hall conductance) 
oi G = /h, different from the well-known 



G = u— (5.12) 
h 

valid for a fractional quantum Hall state (z/ = 1/m is the filling factor). To repair this problem 
one makes the hypothesis that instead of eqs.( [4.15| ) and ( |4.2CI| ) one has 



p(x) = _}l^^p_ ^ ^-(a;) = Mv^H+(x) , (5.13) 

TT OX 

where the subscripts indicate projection on right-going states. With these definitions following the 
calculations of the previous subsection one now straightforwardly reproduces the correct result, 
eq.( p.l2| ). The appearance of the factors y/D in eq.( p.l3| ) indicates that the objects occupying the 



states in figure ^ are not free electrons but rather strongly affected by the physics of the bulk of 
the samples. A more detailed derivation, starting from a Chern-Simons field theory for the bulk 



physics, has also been given by Wen |70 



Beyond reproducing the correct value of the Hall conductance, the above hypothesis leads to a 
number of interesting conclusions. Consider first the creation operator for a real electron (charge 
e) on the edge. Following the arguments of section |4.1|, because of eq. (^.13|) , the bosonized version 



of the electron operator now must create a jump of of height vr/^^, rather than of height tt. 
This leads to 

i)+{x) ^ e-'"t>+^^y^ , (5.14) 

X being the coordinate along the perimeter of the sample. Now, these operators obey the relation 

%l)+{x')%l)+{x) = e^^"/>+(x)7/'+(x') . (5.15) 

But the real electron is still a fermion, i.e.'0_|_ must obey anticommutation relations. Thus m = l/u 
has to be an odd integer. One thus reproduces one of the fundamental facts of the fractional 
quantum Hall effect. From eq. ( |5.14| ) one also finds a decay of the single-electron Green function 



as 

Another fundamental property of the quantum Hall state appears when one considers the 
fractionally charged elementary excitation of charge eu at the edge. A charge-ez/ object is created 
by 

V^+,(x) ^ e-'^^+(") , (5.17) 
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leading to a slow decay of the corresponding Green function, with exponent z/, instead of 1/z/ in 
eq.( p.l6D . One now has the relation 

V^+,(a;')^+.(a;) = e±-'^^+,(a;)V'+.(a;'), (5.18) 

i.e. exchanging the fractionally charged objects one obtains nontrivial (7^ ±1) phase factors and 
these quasiparticles thus also obey fractional statistics. 

A single hypothesis, the insertion of the factors in eq. (|5.13|) , thus reproduces two of the 
fundamental facts about the fractional quantum Hall effect! In addition one obtains a basis 



for the treatment of transport phenomena mentioned at the end of sec. ^.2.1| and results for the 
asymptotics of Green functions. 

5.2 Disorder 

5.2.1 Effects of isolated impurities 

The infinite conductivity in the ideally pure systems considered up to here is a natural but hardly 
realistic result: any realistic system will contain some form of inhomogeneity. This in general leads 
to a finite conductivity, and in one dimension one can anticipate even more dramatic effects: in a 
noninteracting system any form of disorder leads to localization of the single-particle eigenstates 



[[72| , [7^] . How this phenomenon occurs in interacting systems will be discussed in this and the 



following section. 

Following Kane and Fisher |^ , consider first the case of a single inhomogeneity in an otherwise 
perfect one-dimensional system. The extra term in the Hamiltonian introduced by a localized 
potential v{x) is (for spinless fermions) 

-f^barricr OC J (ix^"^ (x) ^'(x) . (5.19) 

Decomposing the product of fermion operators into right- and left-going parts, one has 

7/)t^ = ^p^_^^_^_ + ijji^p_ + + . (5.20) 



In the bosonic representation, the first two terms are proportional to dxcf) (see eq.( 4.15|) ), and 



therefore the corresponding contribution in eq.( p.l9| ) can in fact be eliminated by a simple unitary 



transformation of (p. These terms represent scattering with momentum transfer q -C 2fcF, i.e. 
they do not transfer particles between k-p and —k-p and therefore do not affect the conductance in 
any noticeable way. On the other hand, the last two terms in eq. ( ^.20|) represent scattering with 



|g| 2kp, i.e. from the + to the — branch and vice versa. These terms certainly are expected 
to affect the conductance, because they change the direction of propagation of the particles. The 
bosonic representation of these terms is 

^^barrier = ^^^^ COS 20(0) , (5.21) 

where the potential V{x) is assumed to be centered at x = 0. For this reason, only the value of 
the at X = intervenes in eq. ( p. 21] ). 



One now can integrate out all the degrees of freedom away from x = 0, to obtain an effective 
action implying only the time-dependence of 0(0). Then a renormalization group equation for 
V = V{2kp) can be found as 

^ = (l-A-)^. (5.22) 
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where E = Eqb ^, Eq is the original cutoff, and E is the renormahzed cutoff. 
From eg. ( ^.221 ) it is clear that there are three regimes: 

1. For K > 1 one has V{£ — ^ oo) = 0, i.e. as far as the low-energy physics is concerned, the 
system behaves like one without the barrier. In particular, the low-temperature conductance 
takes the "pure" value G = e^K/h, with corrections of order -y^^ note that in this 
case superconducting fluctuations dominate, and the prefect transmission through the barrier 
can be taken as a manifestation of superconductivity in the one-dimensional system. 

2. For the noninteracting case K = 1, V is invariant, and one thus has partial transmission 
and a non-universal conductance depending on V. 

3. For K < 1, V{i) scales to infinity. Though the perturbative calculation does not provide 
any direct way to treat this case, it is physically clear that the transmission and therefore 
the conductance should vanish. 

Note that the non-interacting case is marginal, separating the regions of perfect and zero trans- 
mission. These results are very similar to earlier ones by Mattis and by Luther and Peschel 
|176| who treat disorder in lowest-order perturbation theory. 

The case of < 1 can be further analyzed considering the case of two finite Luttinger liquids 
coupled by a weak tunneling barrier, as would be appropriate for a strong local potential. The 
barrier Hamiltonian then is 

i^barrier = t (^'1(0)^2(0) + 4(0)^i(0)) ^ — COS 2^(0) . (5.23) 

Here il'i,2 ^-'^^ the field operators to the left and to the right of the barrier. The operators have to 
satisfy the fixed boundary condition ipi (x = 0) = 0, different from the periodic boundary conditions 
we have used so far. Noting that the ipi can be decomposed into left- and right-going parts as 
'ipi = "ipi^ + i'i-, this can be achieved, using eq. (^4.16|) , by imposing the fixed boundary condition 
= 0) = 7r/2 on the boson field fT^, |78[. 
One can now proceed in complete analogy to the weak-l^ case to obtain the renormalization 
group equation 

Again, there are three different regimes: (i) for K > 1 now t{£ oo) — >■ oo, i.e. the tunneling 
amplitude becomes very big. This can be interpreted as indicating perfect transmission, e.g. 
G = e'^K/h] (ii) the case K = 1 remains marginal, leading to a t-dependent conductance; (iii) 
for K < 1 t scales to zero, there thus is no transmission, and G = 0. The results obtained in 
the two limiting cases of small V (weak scattering) and of small t (weak tunneling) are clearly 
compatible: e.g. for K < 1,V becomes large, i.e. at sufficiently low energies one expects essentially 
a tunneling type behavior, and then from eq. (|5.24|) the tunneling amplitude actually does scale to 



zero, giving zero conductance in the low-energy (or low-temperature) limit. For K > 1 a similar 
compatibility of the two limiting cases is found. The global behavior can be represented by the 



"phase diagram" in figure |T0[ For electrons with spin but spin-independent interactions, results 
are very similar: the separation between zero and perfect transmission is at Kp = 1, with Kp = 1 
again the marginal case. In the transmitting region the conductance is G = 2Kpe^ /h. 



These considerations can be generalized to the case of two barriers |74, u%. In particular 



assuming that there are two identical, weakly scattering barrier at ±(i, the effective scattering 
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V 



G = 



J L 



G = Ke ?h 



K 



Figure 10: "Phase diagram" of a localized inhomogeneity in a spinless Luttinger liquid, characterized 
by an exponent parameter according to [^]. The scaling trajectories calculated for weak y or t are 
indicated by arrows. It is clearly plausible to assume direct scaling from weak to strong coupling in the 
whole range of K. 



potential becomes Kff(?) = 2v{q) cos{qd/2). Though in general this is non-zero when V{q) is 
non-zero, for particular values of kp, so that cos(fcF'^) = this potential vanishes, giving rise to 
perfect transmission even for K < 1. This resonant scattering condition corresponds to an average 
particle number between the two barriers of the form u + 1/2, with integer z/, i.e. the "island" 
between the two barriers is in a degenerate state. If interactions between the electrons in the 



island are included, one can recover the physics of the Coulomb blockade [^4], [79 



For the chiral Luttinger liquid, discussed in section |5.1.3| , backscattering events a priori seem 
to be excluded because all the particles are moving in the same direction. In that sense the 
chiral Luttinger liquid can be considered as "perfect" . However, if the quantum hall device has a 
constriction that brings the two edges close to each other, scattering from one edge to the other 
becomes possible and is the equivalent of backscattering. Then similar considerations as made for 
the single-impurity case are possible [0, and in particular the crossover function describing the 
conductance through a resonance as a function of temperature [BTI, B^f and the liV) characteristic 
331 have been obtained. 



5.2.2 Anderson localization of one— dimensional interacting fermions 

The discussion of the previous section was concerned with the effect of at most two impurities, 
weak or strong. Clearly, in that case the effects of coherent scattering from many impurities, which 
typically give rise to Anderson localization, are absent. We now turn to this more complicated 
case which had been studied in fact well before the single impurity work |Q, |S^, . 

In the absence of electron-electron interactions, localization effects can be discussed in the 
framework of a scaling theory |]73[ . Under the assumption that at some short length scale one has 
elastic scattering of electrons off impurities, this theory leads to the following /^-function for the 
variation of the conductance with linear dimension L: 



PiG) 



dln{G) 
dln{L) 



d 



2-- + ... 



(5.25) 



where a is a constant and d the spatial dimensionality. In particular in one dimension this leads 
to a conductance decaying exponentially with the length of the system, exhibiting clearly the 
localized character of the all single-electron states (a fact first shown by Mott [|72| and studied in 



great detail since [^, p9|). 
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In an interacting one-dimensional system (as described by the Luttinger liquid picture of 
the previous section) now a number of questions arise: what is the influence of disorder on the 
phase diagram obtained previously? What are the transport properties? Can one have true 
superconductivity in one dimension, i.e. infinite conductivity in a disordered system? To answer 
these question we discuss below the generalizations necessary to include disorder in our previous 
picture. 

We start by the standard term in the Hamiltonian describing the coupling of a random potential 
to the electron density 



j dx V{x - Ri)p{x) , (5.26) 



where the Ri are the random positions of impurity atoms, each acting with a potential V on the 
electrons. In one dimension one can distinguish two types of processes: (i) forward scattering, 
where the scattered particle remains in the vicinity of its Fermi point. As in the single-impurity 
case, this leads to a term proportional to d(f)p and can be absorbed by a simple redefinition of the 
(f)p field. The physical effects are minor, and in particular the DC conductivity remains infinite, 
(ii) backward scattering where an electron is scattered from k-p to —k-p or vice versa. For small 
impurity density this can be represented by a complex field ^ with Gaussian distribution of width 
= mViq = 2kYf: 

Ht = Y.[dx (e(a;)^L(^)^L.(^) + h.c.) (5.27) 

This term has dramatic effects and in particular leads to Anderson localization in the noninter- 
acting case . 



From a perturbative expansion in the disorder one now obtains a set of coupled renormalization 
group equations |8j 



dKrr 1 



2 

dy 



2{l -K„)y-V 



dt 

— = [^-Kp-K^-y)V (5.28) 
where V = 2D^a/(nu'^) (ua/up)^'' is the dimensionless disorder, y = gi±/{7rUa-) is the dimension- 



less backscattering amplitude, and the are defined in eq. (f4.29|) . These equations are valid for 
arbitrary Ky (the usual strength of bosonization), but to lowest order in V and y. 

As a first application of eqs. ( |5.28| ) one can determine the effect of the random potential on the 
"phase diagram", as represented in figure ||. In fact, there are three different regimes: 

1. for Kp > 2 and gi± sufficiently positive the fixed point is 'D*,y* = 0,K* > 2. Because the 
effective random potential vanishes this is a delocalized region, characterized as in the pure 
case by the absence of a gap in the spin excitations and dominant TS fluctuations. 

2. For Kp > 3 and gi± small or negative one has V* = 0, y ^ — oo, K* > 3. Again, this is a 
delocalized region, but now because y —oo there is a spin gap and one has predominant 
SS fluctuations. 
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3. In all other cases one has V oo,y — oo. This corresponds to a localized regime. For 
small Kp the bosonized Hamiltonian in this regime is that of a charge density wave in a 
weak random potential with small quantum fluctuations parameterized by Kp. This region 
can therefore be identifled as a weakly pinned CDW, also called a "charge density glass" 
(CDG). The transition from the CDG to the SS region then can be seen as depinning of the 
CDW by quantum fluctuations. 

One should notice that the CDG is a nonmagnetic spin singlet, representing approximately a 
situation where localized single-particle states are doubly occupied. Though this is acceptable 
for attractive or possibly weakly repulsive interactions, for strong short-range repulsion single 
occupancy of localized states seems to be more likely. One then has a spin in each localized state, 
giving probably rise to a localized antiferromagnet with random exchange (RAF). 

A detailed theory of the relative stability of the two states is currently missing and would 
certainly at least require higher-order perturbative treatment. The boundaries of the different 
regimes can be determined in many cases from eqs.(|5.2^), and the resulting phase diagram is 



shown in flgure 11. 




Figure 11: Phase diagram of a Luttinger liquid in the presence of a weak random potential (T> = 0.05). 
The full lines represent results obtained directly from the scaling equations ( 5.2^ ), the dashed lines are 



qualitative interpolations. The dotted lines are the phase boundaries in the limit D ^ 0. 

The localization length for small disorder can be obtained from standard scaling arguments: 
suppose that a system with some flxed disorder has a localization length ^q- Then in the 
general case one has ^(I^) = ^oe^^^°'^\ where i{T)Q,T)) is the time it takes for the "bare" disorder 
V to scale up to Vq. From this reasoning one flnds, for the case without a spin gap {gi > 0) and 
weak disorder 

e(I?) oc (l/I?)^/('-^-) . (5.29) 

Note that for Kp > 1, i.e. superconducting fluctuations predominating in the pure case ^ is greater 
than the mean free path A oc l/'D, there is a kind of diffusive regime, contrary to the noninteracting 
case. On the other hand, for Kp < 1 one has ^ < A. In the vicinity of the TS-RAF boundary one 
has 
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The analogous results for the case with a spin gap {gi < 0) are 



^{V) oc (5.31) 
e(P) oc expf , I (5.32) 

There are two points to be noted about this result: (i) for Kp = one has ^ oc P"^/^, which is 
the same result as that found for the pinning length of a classical CDW |]9T|. (ii) the results ( p. 29 



5.30|) and ( ^.31| , ^.32|) are qualitatively different, both in the vicinity of the phase boundaries and 



in the localized states. The transitions are thus in different universality classes, and this strongly 
supports the idea that the localized phases reached through the transition are themselves different 
(RAF or CDG). 

The temperature dependence of the DC conductivity can be obtained noting that at finite 
temperature there are no coherent effects on length scales larger than Vp/T. One therefore stops 
renormalization at e^ = vp/laT). As long as T> remains weak one can still use the Born approx- 
imation to obtain 

= , (5.33) 

where ctq = e'^Vp/27chD^ is the lowest order conductivity. In the delocalized phases one then finds 
a conductivity diverging as cr(T) ~ where 7 = K* — 2 in the TS case and 7 = ii'* — 3 in the 

SS case. On the phase boundaries one has universally a ~ 1/T. In the localized region V diverges 
at low temperatures, and a perturbative calculation thus becomes meaningless. However, the 
conducting-localized crossover can still be studied at not too low temperatures In particular. 



the high-temperature conductivity is found to vary as a ~ ^f. This is the perturbative result 
first found by Mattis [|7S[ and by Luther and Peschel and also reproduced by the single- 



impurity calculations |7^. The high-temperature behavior thus can be understood in terms of 
scattering off the individual impurities. On the other hand, at lower temperatures one necessarily 
comes into the region where V increases sharply. This has its origin in coherent scattering from 
many impurities and ultimately gives rise to localization. 

One can finally notice the effects of different types of interactions on localization. Roughly 
speaking, for forward scattering repulsion ((72 > 0) enhances localization whereas attraction weak- 
ens it. In particular, strong attraction leads to vanishing effective random potential. The delocal- 
ized state then can be considered to be a true superconductor in the sense that there is infinite 
conductivity even in an impure system. The effect of backward scattering interactions is opposite 
to that of forward interactions. 



5.3 The spin— 1/2 chain as a Luttinger Uquid 

One of the fundamental models of solid state physics is the Heisenberg model of insulating magnets. 
In the one-dimensional case ( "spin chains" ) its Hamiltonian takes the form 

1=1 

L .1 



1=1 
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Here Si = (S'f , S'f , S'f ) is a spin operator (we will first concentrate on the case of spin 1/2 so 
that S; • S/ = 3/4) acting on site /, A is an anisotropy parameter that allows one to treat the 
antiferromagnetic (A = 1), the ferromagnetic (A = —1), and general anisotropic cases, and 
periodic boundary conditions imply S^+i = Si. We notice that the Hamiltonian conserves the 
z-component of the total spin (for |A| = 1 total spin is also conserved). 

The spin model can be transformed into a model of spinless fermions, noting that and 
Sf anticommute. The Jordan-Wigner transformation then relates spin to fermion operators 
{ai,al) via 

5;+ = a] exp l^ivr ^ a]a^ , S'f = a\ai - i . (5.35) 

Presence or absence of a fermion now represent an up or down spin, and the exponential factor 
insures that spin operators on different sites commute, whereas fermionic operators of course 
anticommute. The transformation can now be used to rewrite the spin Hamiltonian (|5.34| ) in 
terms of fermions as 



1=1 



The "spin-flip" terms thus give rise to motion of the fermions, whereas the S^-S^ interaction 
leads to a fermion-fermion interaction between adjacent sites. 

It is instructive to see in detail how one can pass to the continuum limit and then to a bosonic 
model starting from eq.( ^.36 ). We start by passing to momentum space in the standard way: 



ai = 4fT.(^k^''' , (5.37) 

with the momentum sum restricted to the first Brillouin zone: — vr < k < n. Insertion into 
eg. ( ^.361) then straightforwardly leads to 



H = J2i-cosk-A)aia, + jYl^if'^ + k2-h-h)e'^''-'*^al^^^^ (5.38) 
= cos - A)alaf^ + ^ ^ ^ ~ ^^^^^^ - h + q)) al^+gal^-ga^^ak^ , (5-39) 

k k\k2q 

where the 5-symbol insures momentum conservation modulo 27r: 5{x) = 1 if x = mod 27r, and 
5{x) = otherwise. In order to be close to the standard case of a band minimum at /c = we have 
shifted the origin of /c-space by vr. This amounts to a rotation of every other spin in ( |5.34| ) by 
71 around the z axis. In eq. (|5.39|) the form factor appearing in the interaction has been properly 
symmetrized. This will be useful subsequently. 

At least for weak interaction (small A) only interactions involving states close to the Fermi 
energy are important, and one therefore can then map eq.( p.38D onto the spinless Luttinger model 



discussed in section ^. For the g2 interaction, processes with either ki ^ ^4 or ki ^ k^ contribute 
(the ~ sign is meant to indicate that both momenta are close to the same Fermi point). To cast 
the second type of processes into the form of (^4.8|) one has to commute two fermion operators, 
giving rise to an extra minus sign. The coupling constant at small momentum transfer then is 

g2 = 2A(1 -cos2A;f) . (5.40) 
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Note that this vanishes when the spin chain is nearly fully spin-polarized {kp —>■ 0,7t). 

The determination of the appropriate parameter is slightly less straightforward. Naively, 
one would expect that processes where all four states are in the vicinity of the same Fermi point 
contribute, giving (74 = 2A. However, as pointed out by Fowler [^, this is not correct: in fact, in 
the lattice model there are corrections to the bare fermion energy coming from the exchange part 
of the first-order Hartree-Fock selfenergy, given by 

2A 2A 

Shf(^) = — -j-^io-lo-p) = sin/cpcos/c . (5.41) 

This leads to a renormalization of the Fermi velocity 

2A 

f F = sin kp — > sin kp H sin^ kp . (5-42) 

71 

On the other hand, in the continuum model of section |^, such Hartree-Fock selfenergy corrections 
must not be considered (formally by appropriately normal ordering the interaction terms), oth- 
erwise they would be infinite due to sums over infinitely many occupied states. The proper way 
to account for the finite velocity renormalization, eq. (|5.42|) , then is to include this via a properly 



chosen (^4. In the present case, in order to reproduce (|5.42|) one then has to set 



94 = 92 = 2A(1 - cos 2kp) . (5.43) 

Note that at half-filling this is twice the naive expectation. 

Following the steps of section 11, the low-energy excitations of the spin chain model then are 



described by the Hamiltonian ||94 

H = ldx{ ^H(x)2 + ^{d.<Pr 1 ■ (5.44) 



This is of course our standard Hamiltonian, i.e. spin chains are Luttinger liquids Specializing 
to the case of zero applied field, so that the total magnetization vanishes and one has kp = 7r/2, 
the parameters are given perturbatively by 




u = Jl + — , K = ^== (5.45) 



1 + 4A 



These results agree to first order in A with the exact ones, to be discussed in section |5.3.1| below, 
and are close to them over much of the parameter range. 

For an exactly half-filled band extra umklapp scattering processes are possible, with fci 2 ~ 7r/2, 
^3,4 ~ ^7r/2, and vice versa 0, 0. In a continuum representation, these operators become 



H^ = Ajdx ((4-9.4) {^Ld.^L) + h.c.) = J cos40(x) (5.46) 

This operator has scaling dimension 4K and therefore is strongly irrelevant for small A. However, 
it becomes relevant for < 1/2 and then in particular is responsible for the creation of a gap 
in the excitation spectrum and long-range antiferromagnetic order in the case of an Ising type 
anisotropy (A > 1). 
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5.3.1 Physical properties of the spin 1/2 chain (small A) 

From the bosonized form of the Hamiltonian one immediately obtains the specific heat of the spin 
chain as 

The susceptibility for a field applied along the z-direction can also be obtained, noting that dxcf) 
is proportional to the fermion density, i.e. the 2;-component of the magnetization: 

X = (5.48) 

TT U 

Beyond these thermodynamic properties the fermionic analogy allows one to study correlation 
functions of the spin chain using the bosonization formalism developed in section |4.1| . We start 
by rewriting the fermion operators in a continuum form: 

ai = tpnix) + tfjiix) , (5.49) 

where x = la and in most of what follows the lattice constant a is set to unity. In order to 
calculate spin-spin correlation functions we further need a representation of the exponential factor 
in eq.( 5.35| ). This is given by 

exp l^ivr^^^a^a^.j =e'^^^e-''^(^) , (5.50) 

where the first exponential factor comes from the mean value of the density and the second one 
represents the effect of fluctuations about that mean, see eq.( |4.15| ). 

We can now use the continuum representations ( |5.49| ) and ( |5.50| ) together with the fermionic 
representation of the spin-1/2 operators, eq.( |5.35D , to give a bosonic representation of the spins: 

St = ^'I^4^r{x) + ipt^L{x) +^t^R{x) + tp^^L{x) = --d^(j) + — cos{2(l){x) -2kFx) 

TT TTa 



7ra 



Here we have omitted the Klein factors which will be unimportant in the following, and have kept 
a finite a as a short-distance cutoff, to be taken of the order of a lattice constant. Further, in 
we have restored the factors (—1)' that were lost in going from eq. (|5.36|) to eq. (|5.39|) . Using these 
expressions we can now follow the calculations of section [4.1| to obtain the spin-spin correlation 
functions. In particular, at zero temperature |94] 

IT <iH^\<iH^\\ ^ x^-u'^T^ Acos{2kFx) 



(T.^+(r)5o"(0)) = ^1 cos((2A:f - 7v)x) ^''^2^2^kIi/4K+i + ^^2 + X^^^"^ ' ^^'^^^ 

Here the constants A, Bi cannot be reliably determined by the present methods, however in 
the nonoscillating part of the z-z correlation function there is no undetermined parameter, and 
very recently B2 has been determined. |^8| In most cases the dominant contributions in eq.( |5.52| ) 
come from the second, quickly oscillating terms. The alternation indicates the expected tendency 
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Figure 12: Real (a) and imaginary (b) part of the quantum critical scaling function 
isotropic Heisenberg antiferromagnet {K = 1/2). 



(x, y) for the 



towards antiferromagnetic order, but correlations do decay (with a rather slow power law), and 
there is thus no long-range order, as to be expected in one dimension. 

Fourier transforming eq. ( p.52| ) gives, after analytic continuation to real frequencies and gener- 
alizing to nonzero temperature [B^ 
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at T = 



(5.53) 



Here B{x,y) = T{x)T{y)/T{x + y) is Euler's beta function, q' = q — ^kp, and the last term has 
been introduced "by hand" in order to reproduce the known logarithmic result for K = 1. For 
the transverse (H — ) correlations, an analogous result holds, with K 1/{4K) and q' = q — tt. 
One should notice that eq.( p.53| ) is valid uniformly for arbitrary ratios between the energies T, 
u, and u\q'\, as long as all of them are small compared to the exchange energy (which here play 
the role of the bandwidth). At zero temperature, a finite imaginary part of x\\ and x± exists in 
the V-shaped region \uj\ > u\q'\. This represents precisely the regions where the Bethe ansatz 
solution produces the two-spinon continuum. | |100| | The magnetic scattering cross section obtained 
from eq.( |5.53|) is in excellent agreement with experiment. ||101| , |102| , |103 | 



At finite temperature, the susceptibilities take the general quantum critical scaling form [|104 

A, , , 



UJ 



uq 

4^' 4^, 



(5.54) 



given by the products of beta functions in eq.( ^.53|) . Figure [12 



J^2-2K^\\ 

with the scaling functions 
shows plots of this function for the particularly relevant case of the isotropic antiferromagnet, 
K = 1/2. Note that these function has maxima for x ~ For x ~ ?/ — > cxd these maxima develop 
into the (zero-temperature) square root singularity of eq. (|5.53|) . 
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Generalizing the arguments leading to eq. (|5.51|) it has been shown |[L05|| that the nearest- 
neighbor interaction contains an oscillating part of the form 

Si ■ Si+i ^ — sin(20(x) - 2kFx) . (5.55) 
na 

Correlation functions of this operator, which describes the tendency toward dimerization (or spin- 
Peierls order), then decay with the same power law as those of the z component of the spin. 

An important consequence of the above discussion is that the low temperature correlations 
are only determined by two parameters: one, K, determines the power laws of the decay of 
correlations, the other one, u, is the velocity of the excitations. In the previous section, we have 
given explicit formulae for these constants, eq.( |5.45| ), however these can only be expected to be 



valid for small A where the linearization of the fermion spectrum is expected to be quantitatively 
reliable. It is clearly interesting to determine these parameters outside the perturbative regime. 
The exact Bethe wavefunction is, at least up to now, too complicated as to allow calculation of 
correlation functions (however, for some recent progress see ||106| , |107|| ). One point where a non- 



perturbative result can be obtained is the isotropic Heisenberg antiferromagnet A = 1 in zero 
field, so that kp = 7r/2. Then spin rotation invariance requires the longitudinal and transverse 
correlation functions to be equal, and thus i^" = 1/2 is needed in eq. (|5.52| ). Note that this implies 
that the dimerization correlations described by eq.( |5.55D also decay with a 1/r power law. 



For general A a more indirect approach to the determination can be used [^. One notices 
that the bosonized Hamiltonian leads to a variation of the ground state energy with the number 
of right- and left-going particles as 

E{N,J) = — [uKj' + -N'j , (5.56) 

where N and J are the sum and difference of the number of added right- and left-going particles. 



and eqs.( |4.10| ) to ( [4.12 ) have been used. Haldane introduces the charge and current velocities 



vn = u/K and vj = uK. On the other hand, the variation of the ground state energy with 
and J can also be obtained from the Bethe solution (or even from a numerical solution on a finite 
lattice, if an exact solution is not available). In particular, for the spin-1/2 chain in zero field 



E(N, 0) = ('^-^rccosAj.^i^iV; ^ ^^^^^^ 

2 arccos A L 

This fixes u/K. Note that the coefficient of A^^ in eq.( |5.56| ) is proportional to the inverse sus- 
ceptibility. The variation with J can also be obtained, noting that J oc d^O can be related to 
the stiffness constant of the x and y components of the spin (compare eq.( |5.51| )). However, it is 



simpler to use directly the known result for the velocity of the elementary excitations ||109|| : 



vrVl - A2 

u = — . 5.5J 

2 arccos A 



Alternatively, the velocity can also be obtained from the low-temperature specific heat [p. 10 
leading to the same result. This then fixes the correlation exponent as 

K=- — , (5.59) 

2[7i — arccos A) 

first obtained in |Q by a scaling argument using the exact solution of the fully anisotropic XYZ 



model [|111|| . The exact results (|5.58|) and (|5.59|) agree to linear order in A with eq. (|5.45|) , as 
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expected, and even beyond first order the approximate results are quite close to the exact ones. 
One should also notice that to first order in the interaction A the coefficients in eg. ( p. 561 ) can be 
calculated using the zeroth-order wavefunction, i.e. the noninteracting Slater determinant. This of 
course reproduces eq. ( |5.43| ) and probably represents the "cheapest" way to calculate the anomalous 
exponent K. The exact expressions ( p.58| ) and ( p.59| ) are valid for vanishing magnetization, i.e. no 



applied field. Numerical results for a finite magnetization along z (i.e. a non-half-filled fermion 



band) have been given by Haldane pS 



5.3.2 The isotropic antiferromagnet (A = 1) 

Up to now we have ignored the effect of the umklapp operator, eq. (|5.46|) . Indeed, as long as 
it is irrelevant, i.e. in the region of planar anisotropy A < 1, umklapp interactions only lead to 
subleading corrections to correlation functions ||38|. However, at the isotropic antiferromagnetic 
point A = 1, the umklapp operator is only marginally irrelevant. Then the corresponding coupling 
constant is renormalized as ||112| , |38| , |40[| 



% - -^W . (5.60) 
with solution 

^ \ + T[hg\n.L IuL/Lq ^ ^ 

Here h = 4/^3 is a normalization constant fixed by the requirement that the correlations of the 
marginal operator decay as with unit coefficient, L = e^ is the rescaled short distance cutoff 
(lattice constant), and Lq = exp{—l/{7ibg)). The determination of the bare coupling constant 
g is not entirely trivial. Naively from eq.( |5.46D one would set it of order A. However, for the 
isotropic antiferromagnet A = 1, giving a coupling of order unity which would be expected to be 
outside the perturbative domain of validity of eq. ( p.60| ). A precise determination can be achieved 
noting that the marginal operator affects the low-lying excited states in a finite system in a well- 
understood way involving only the combination Tthg{t) ||112| , One finds in particular that the 



predicted linear variation of 7rbg{i) with ln(L) is satisfied to within a few percent, and an order 



of magnitude better agreement is found if two-loop corrections are included | 113 |. From the two- 
loop calculation the bare coupling constant (at L = 1) can be estimated as g ^ 0.24. A priori, this 
would be expected to be an effective coupling, reproducing correctly the long-distance behavior, 
however in fact the spectra of all but the shortest chains are rather well fitted by this value. On 
the other hand, there are considerable uncertainties associated with this estimate: fitting the same 
data to the one-loop function one finds g ~ 0.11. 

The umklapp operator lead to logarithmic corrections in the temperature dependence of the 



spin susceptibility, of the form ||114 | 



'^(^) = ^ + 21^) • 

with To ~ 7.7, obtained from fitting eq. ( |5.62| ) to exact finite-temperature Bethe ansatz results. 
Experimental results on the nearly perfectly one-dimensional compound Sr2Cu03 [[L15|| are in 
excellent agreement with the logarithmic law ( p.62| ) and with Bethe ansatz results ||1 14|| over 



a wide temperature regime. Similar logarithmic corrections also exist in the zero temperature 
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magnetization curve M{H) ||108|| . On the other hand, there are no such corrections in the low- 
temperature specific heat. 

The marginal operator further produces multiplicative logarithmic corrections in various cor- 
relation functions, e.g. [^, ^ ^ 



(T.S,(r) • So(0)) 



Ln 



(5.63) 



where the prefactor has only been determined very recently. [|116|] An identical logarithmic correc- 
tion is found in correlations of the dimerization operator, eg. ( p. 551 ), however with exponent —3/2 
instead of 1/2, i.e. the dimerization fiuctuations are logarithmically suppressed compared to the 
antiferromagnetic ones. Also the numerical prefactor in those correlations is not currently known. 
Similarly, the staggered susceptibility has a logarithmic correction factor: 



A 



x{n,0;T) = -\n'/\To/T) 



T 



(5.64) 



Numerical investigations of correlation functions initially shed doubt on the existence of these 
logarithmic corrections [ |1 1 7|| , however more recent numerical ||118| , |119| , |120| , |121| , |122|| and analytic 
106| , |107|| work provides ample evidence for their existence. In particular, the constants in eq.( |5.64| ) 



have been determined as A = 0.32 ± 0.01, Tq = 5.9 ± 0.2 [ 122|| . The difference of Tq with the value 
obtained from the susceptibility is possibly due to the fact that the susceptibility data were fitted 
in the asymptotic low-temperature region (T < 10~^), whereas the staggered susceptibility was 
calculated at somewhat higher temperatures. 

Apart from being directly accessible in neutron scattering experiment s , [p. 1| , |102| , |103|| the 
above staggered spin correlation function leads to characteristic temperature dependences in the 
longitudinal and transverse NMR relaxation times ||121|| (under the assumption of a hyperfine- 
dominated relaxation) 



oc ( In 



T 



1/2 



T, 



oc 



2G 



T T, 



(5.65) 



in good agreement with experiments on Sr2Cu03. ||123|| In the absence of the marginal operator 



the logarithmic factors would be absent ||124|| , leading in particular to a temperature-independent 
Ti. 

We finally mention that for A > 1 the spins are preferentially aligned along the 2;-direction, 
and one then has a long-range ordered ground state of the Ising type. There is thus a phase 
transition exactly at the isotropic point A = 1. In the fermionic language, this corresponds to a 
metal-insulator transition. ||125|| 

We can summarize this section by noting that quantum spin chains provide one of the experi- 
mentally best established cases of Luttinger liquid behavior. This largely due to two facts: (i) the 
relative ease with which one can define a microscopic Hamiltonian for a given experimental sys- 
tem: one rarely has to go beyond a slightly modified Heisenberg model, with very few interaction 
constants to be determined; this is to be compared with the difficulties one encounters in conduct- 
ing systems: long-range interactions, electron-phonon interactions, etc.; (ii) the availability of a 
large number of experimental techniques which give results directly comparable with experiment 
and, concomitantly, the possibility of using either very well-controlled (often exact) theoretical 
methods or numerical approaches which for spin systems are much more reliable than for itinerant 
fermions. 
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6 Spin ladders and coupled Luttinger liquids 



6.1 Coupled spin chains 



Investigating models of coupled parallel chains is of interest for a number of reasons: (i) quasi- 
one-dimensional antiferromagnets always have some form of interchain coupling, usually leading 
to three-dimensional ordering at sufficiently low temperatures; ||126| , |127| , |128|| (ii) there is a number 



of "spin-ladder" compounds containing a small number of coupled chains; [|129| , |130|| (iii) coupled 



spin-1/2 chain models can be used to describe higher spin quantum numbers. 1 131 ] 

Consider N coupled spin-1/2 chains with spin degrees of freedom S^, j = 1..N, described by 
the Hamiltonian 



N 



H = Y^ H{Sj) + XjkHc{Sj, Sfc) 



HriSj, Sfc) 



k,i 



(6.1) 



j<k 



where i labels sites along the chains, j,k label the chains, and H{Sj) is of the form ( [5.341 ). For 
Xjk = —1 the ground state of this model is exactly that of the spin-A^/2 chain, each site being in 
a state of total spin A^/2, ||13i[| but the model can be considered for general A, both ferromagnetic 
(A < 0) and antiferromagnetic (A > 0). Performing now the Jordan-Wigner transformation for 



the S, separately and going to the boson representation the Hamiltonian becomes] 132 



H 

+ 



dx 



TTUK ^ 



-U(xY + 



^ ^ /" rfx{AijfcCos(2(0j + 0fc)) + A2,jfccos(2(0j - 0^)) + A3jfcCos(^j - Ok)} , (6.2) 
y'^'^) j<k-' 



Here 



= (01, 02, ••, ^Tv), = J2j (pj/ yN, the constants u,K,K all depend on the different con- 
stants in the original Hamiltonian, and the Ai,2,3;jfc are all proportional to Xjk- 

Elementary power counting, using the result eq. (|5.63| ) for the spin correlations, shows that the 
coupling term always is a strongly relevant perturbation. Consequently, an explicit renormal- 
ization group calculation ||132|| shows that either A2 or A3 always scale to strong couphng, i.e. the 
"relative" degrees of freedom 0j — 0^ all acquire a gap. In particular, for not too strong anisotropy, 
the A3 operator dominates, giving rise to long-range order in the 6j — 9k and correspondingly ex- 
ponential decay of the 0j — 0fc correlations. Integrating out these massive degrees of freedom an 
effective Hamiltonian for the "global" mode is found: 




T^KUixY + 



'kK 



(40) +^cos(/iVA^0) 




(6.3) 



where /x = 2 for even and /i = 4 for odd A^, the coefficients m, g are renormalized, and h 
is an external magnetic field applied along the z direction. Similarly, the leading contribution to 
spin correlations comes from the operators 



S+(x 



S\x) oce*"^cos(2ViV0) 



(6.4) 



where the second equation applies to odd only. 

From this a number of important conclusions can be drawn. |p.32|| We first notice that massless 
excitations and the corresponding slow algebraic decay of correlation functions are only possible 
if the cos term in eq.( |6.3|) is irrelevant. Moreover from eq. (|6.4| ) it follows that spin correlation 
functions are isotropic only if i^' = 1/(2A^), implying a decay as 1/x, as in the 5* = 1/2 case. For 
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the case of odd (equivalently, for ferromagnetic A, for half-odd-integer S) this is indeed the 
correct behavior: for /i = 4 the cos term is marginally irrelevant. Thus both antiferromagnetic spin 
chains with half-odd-integer S and for odd numbers of coupled S = 1/2 chains massless behavior 
is predicted, with correlations asymptotically decaying like those of a spin-1/2 chain. [11^] There 
is both numerical ||133| , p.34| , |135| , |136|| and experimental [p.29|| evidence that this is correct. 

On the other hand, for even N (equivalently, integer S) the cos term in eq.( |6.3D is strongly 



relevant and therefore generates a gap in the spin excitations. For the integer-S* spin chains 
this is the well-known and verified Haldane prediction ||137|| , but there also is a gap for any 
even number of coupled chains. [T^, |138|| Analogous conclusions concerning qualitative differences 
between even and odd numbers of coupled chains have also been reached based on the nonlinear 



cr-model. ||139|| The gap implies exponential decay of spin correlations. Numerical [|135| , |136|| and, 
at least for N = 2, experimental work [|129|| confirms this picture. Another prediction, again valid 
both for integer-S* antiferromagnets and even numbers of coupled chains, concerns the effect of an 
applied magnetic field: 1 132 , 140 , 141 1 as long as the field is smaller than a critical field he (x Ag, 
the ground state is unchanged and has zero magnetization. However, beyond he the magnetization 
is expected to increase as M oc \^h — he- Experiments on 5* = 1 antiferromagnetic chains confirm 
this prediction. ||142 1 

A natural question left open by the above considerations concerns quasi-one-dimensional anti- 

p6i WH or Sr2Cu03,|l2|, |T| which can be considered as the A^ ^ oo 



ferromagnets like KCuFg 
limit of the above model. One clearly expects (and observes) true antiferromagnetic order at suf- 
ficiently low temperatures, at first sight in contradiction both with the exponential decay of spin 
correlation predicted for even A^ and the universal 1/x law for odd A^. However, one should note 
that on the one hand the correlation length of the even-A^ systems is expected to increase quickly 
with increasing A^, and that on the other hand the 1/x correlation law of the odd-A^ systems 
also is expected to be only valid beyond a correlation length .^(A^) which increase with A^. In the 
thermodynamic limit N ^ oo this then is perfectly consistent with the existence of long-range 
order. Theoretical treatments of magnetic order in quasi-one-dimensional antiferromagnets can 
be based on a mean-field treatment of the interchain interaction ||143| , |144|| which gives quantitative 
predictions for systems like KCuFa or Sr2Cu03. 



6.2 Two coupled Luttinger liquids 



It is clearly of interest to see what happens to the peculiar one-dimensional behavior when one puts 
chains in parallel. This question is of relevance for the understanding of quasi-one-dimensional 
conductors, doped spin ladders, il46|| few-channel quantum wires, 0] and generally for the 
understanding of possible non-Fermi-liquid behavior and correlation-induced superconductivity 
in higher-dimensional solids. Of particular interest is the effect of an interchain single-particle 
tunneling term of the form 



+ h.c.) 



(6.5) 



where ip^^j is the fermion field operator for right (r = +) or left (r = — ) going particles of spin s on 
chain i. Simple scaling arguments ||147|| lead to the "phase diagram" shown in fig.|l^. The dashed 
line represents the crossover below which single-particle tunneling becomes strongly relevant and 
below which one thus expects Fermi liquid like behavior (an alternative interpretation is due to 
149|| ). The full lines indicate where two-particle or particle-hole hopping becomes 



Anderson ||148| , 

relevant. The most plausible interpretation is that this is the temperature where three-dimensional 
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Figure 13: Qualitative phase diagram in the temperature-ii'p plane for Luttinger liquids coupled 
by interchain hopping. K± = 3 ± -\/8- 



long-range order of some type sets in. A more detailed discussion of this in quasi-one-dimensional 
systems (a thermodynamically large number of chains) has been given elsewhere. [Il3|, |150|| It is 



also worth noting that a different approach to the crossover to higher dimensions, working with 
continuously varying spatial dimension, comes to similar conclusions. [|1 5 1|| 

As a model for doped spin-ladder systems, as well as a first step towards a many-chain system, 
one can study the two-chain case. [ 152 , 153 , 154] The tunneling term, eq.( |6.5|) , then leads to a 



splitting of the single-particle bands into symmetric and antisymmetric combinations which we 
label by transverse wavenumbers k_i_ = 0, vr. Now each kj_ mode can be bosonized separately. 
Introducing the linear combinations (l)^± = {(p^o ± (f)i^T^)/\^ (u = p, a) the Hamiltonian (including 
t_i_) then takes the form 



H 



Hq + Hint,l + H,, 
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Hn 



dx 

' — n.rr ^ 
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2(7ra)2 



(ia;{cos20o-, 
1 



(cos26'p„ + cos20o-_ — cos2^^o-_) — cos26'p_ cos26'o-_} 



(2) 

7 



7=± 



P7 



Here g^^ = fifoooo + 19owtto^ g'^',^^^ is the coupling constant for an interaction scattering two 
particles from -states {a,b) into {d,c). The signs of the different interaction terms in eq. (|6.6|) 
have been determined following the reasoning explained in the appendix. Here we consider a 
case where there is only intrachain interaction, implying that all bare coupling constants are 
independent of k±. 

(2) 

For the pure forward scattering model {gi = 0) the only nonlinear interaction ((^ootttt) scales to 
infinity, leading to a gap in the (p— ) modes and in half of the (cr— ) modes. The remaining (cr— ) 
modes are protected by the duality symmetry under ^ 9o--, the (a—) sector is in fact a critical 
point of the Ising type for gi = 0. |p.54 ] For repulsive interactions here the dominant fluctuations are 



(2) 

2(7ra)2 

(2) 



J dx cos 29p_ (cos 20O— + cos 29o 



(6.6) 
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of CDWtype, with decay proportional to r (3+2ii'p)/4^ and K'^ = {TTVF—g2+gi/2)/{Trvp+g2—gi/2). 
This state can be labeled by the number of massless modes as ClSl^, where quite generally CnSm 



denotes a state with n massless charge and m massless spin modes. ||155 

For nonzero gi all interactions scale to strong coupling, only the total charge mode remains 
massless (CISO), reflecting the translational invariance of the system, and all spin excitations 
have a gap. The physics in this regime can be determined looking for the semiclassical minima 
of the different cos terms in eq.( p.6D . One then finds that the CDW correlations now decay 
exponentially, and for the interesting case gi > 0, (72 > fi'i/S, corresponding to purely repulsive 
interaction, the strongest fluctuations are now of "d-type" superconducting pairing, [T^, |154|| 
with decay as r~'^/'y'^^p\ Labeling this state as "d-type" seems appropriate because the pairing 
amplitudes at = and vr intervene with opposite sign. In real space, this corresponds to pairs 
formed by two fermions on the two different chains. Note that even for weak interactions where 
Kp — i> 1 this decay is very slow. The Akp component of the density correlations also has a power 



law decay, however with an exponent 2i^p, much bigger than the SCd exponent. ||154|, |155|| The 



full phase diagram in the gi-g2 plane is shown in fig.|T^. For gi < Q the diagram is identical to 
the single-chain case, however, for (71 > the behavior is changed dramatically, and in particular 
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Figure 14: Phase diagram of the two-chain model. The different dominating fluctuations are: 
SCd: "d-type" pairing; OAF: orbital antiferromagnetism; SS: "s-type" pairing; CDW: charge 
density wave. The critical lines (71 = and gi = 2g2 are of Ising type. 



superconductivity is predicted for repulsive interactions, for example for the Hubbard model which 
would be represented as gi = g2 in the present language. 

Remarkably, results basically identical to this weak coupling analysis can be obtained assuming 
strong repulsive interactions in the individual chains, so that one is for example in the regime where 
interchain electron-hole pair hopping is more relevant than single particle tunneling, Kp < 1/3 in 



fig.|T3|). It is then more appropriate to bosonize the degrees of freedom of individual chains, rather 
than working in fc^-space. Under renormalization one then generates an interchain interaction of 



the form [^, [156 1 



jk 



X 



(27ra)- 



cos 



{jx[cos(V2(^,, - e^k)) + cos(V20,j) cos(y20,fc)] +Vsm{V2^^j) sm{V2(P„k)]} ,(6.7) 
where j, k now are chain indices. Remarkably, this term leads to properties identical to those 



found in weak coupling. [|154|] First, for gi = one has V = J±. Then Hjk is invariant under the 
duality (paj ^ Gaj, and one has an Ising critical theory. Secondly, for (71 > one finds V ^ J±. 
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Then the fields (j)pi — (f)p2, 4>ai +(f>a2, and 6^1 — 0a2 become long-range ordered and one is in a C15'0 
state. Power law correlations again exist for SCd and Akp density fiuctuations, with the same 
scaling relation between the two exponents as in weak coupling. However, because now Kp < 1/3 
the Akp density fiuctuations actually dominate. The equivalent results in the weak and strong 
coupling regime very strongly suggest that the two-chain model is in the same phase for weak and 
strong repulsion. This point is further supported by considering the "t-J ladder model" for strong 
interchain exchange. ||1 5 7| ] 



Numerical work on the two-chain model is in agreement with the existence of d-type pairing, [ [L58| , 



I^ , |16(J| , TWl , 162 | the evidence for the Akp density fiuctuations is however inconclusive. [ 163(] Con- 



cerning experimental observation, one should notice that the SCd state becomes localized by weak 
disorder. ||164| 

6.3 Summary 

In this section we have discussed a number of results, mainly analytical, on the effect of different 
forms of interchain coupling on the Luttinger liquid behavior of strictly one-dimensional systems. 
As far as spin chains are concerned, the most spectacular result is the "oscillation" between even 
and odd numbers of chains, reminiscent of (and formally related to) the Haldane phenomenon 
in spin-S antiferromagnetic chains. There is both experimental and numerical evidence for this 
behavior, as discussed above. 

For conducting chains, we have only discussed the two-chain case. The most interesting 
conclusion here was that for the Hubbard model (and a rather wide class of its generalizations) 
with purely repulsive interactions a d-wave superconducting state is predicted from weak up to 
rather strong repulsion. This seems to be one of the first cases where there is a reliable theoretical 
argument in favor of superconductivity in the (repulsive) Hubbard model. Only for very strong 
repulsion does a. ikp CDW predominate. These results in principle apply directly to doped spin 
ladders, [[L46|| and to few-channel quantum wires. Concerning doped spin ladders, one of the most 
interesting experimental questions certainly is whether the superconducting state is indeed of d 
type as predicted. 

A number of results exist for larger numbers of chains. ||165| , |166| , [167| , |168|| In the perturbative 



region for weak repulsion generally again a d-wave superconducting state is found, however for 



stronger coupling different phases are found. 1 166 1 
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A When to bosonize in peace 



The fundamental ingredient for the following is the expression eg. ( 4. 16 ) for the single-fermion 
field operators 

1 

7l V "^^'""'"^ v^^« 



'^±A^) = —ff XI «±,a,fce'''^ = hrn "^^f— exp {±.%kFX =F %(l)cj{x) + %Qa{x)\ 



(Al) 
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where a can be the spin or some other internal degree of freedom of the fermions and we have 
replaced the [/-operators by the Majorana ("real") fermion operators T]±^cr introduced to guarantee 
proper anticommutation between the V^'s-IPBI They satisfy the anticommutation relation 



[r]r,r]s]+ = '^Sr,s , (A2) 

where r and s are compound indices containing both the chirality ± and the internal degree of 
freedom a. Eg . (|A2|) implies in particular (//r)^ = 1- The Majorana fermions can be represented 
by standard (Dirac) fermion operators as 77^ = cj + c^. Note that there is just one isolated 
fermionic degree of freedom per branch, and that these degrees of freedom will not appear in the 
bosonized Hamiltonian if properly handled. 

Using eq.(|Al|) and its generalization to cases with spin and other "internal" degrees of freedom 
like perpendicular momentum indices in coupled chain problems, a typical fermion interaction 
term becomes 

ipj^ip^l^ip^ipj = rjaTiisri^ris x (boson operators) = hais-ys x (boson operators) , (A3) 

where the second equality defines ha/BjS- This operator, responsible for taking into account fermion 
anticommutation properly, is nevertheless carefully passed under the rug in the vast majority of 
the literature, thus leaving a purely bosonic Hamiltonian to be considered, as implied in the term 
"bosonization" . I will here investigate under which conditions this is allowed. 

The relevant situation is that all the indices 7^ in hap-ys are different from each other. 
Otherwise the anticommutation rule 

[Vr,Vs]+ = 26r,s^Vr = 1 , (A4) 

allows to simplify hap-ys- I will therefore only consider the general case. First note that 

hlp,s = 1 , (A5) 
hai3.ys thus has eigenvalues ±1. Secondly, 

[hal3-yS, hKXnu]± = , (A6) 

where according to whether an even or odd number of pairs of indices taken from the two sets 
(a,/?, 7 5) and (k, A,/i, z/) are equal the commutator (even case) or anticommutator (odd case) is 
to be used. Finally, permutation of indices leads to sign changes: 

It is now clear that if all the /I's occurring in a given Hamiltonian commute, they can be 
simultaneously diagonalized, which means that it will be possible to replace each of the /I's by 
ibl, leading to a purely bosonic Hamiltonian. This clearly is the case if all the /i's occurring have 
an even number of indices in common. In the opposite case some of the /I's do not commute, 
therefore can not be simultaneously diagonalized and not be eliminated from the Hamiltonian. 
Bosonization then is not possible. 

As a simple example consider the single-chain Luttinger model with spin, sec. [4.2| . the four al- 
lowed values of the discrete indices are 1 = (+,|),2 = (+,|),3 = (— ,|),4 = (— ,i). Consequently, 
only /11234 can occur, and according to the eigenvalue chosen the backward scattering interaction 
takes the form ±(71 cos(-\/80o-)- The choice of eigenvalue of /?.i234 affects however the expressions 
for correlation functions: for example /11234 = ±1 implies rjiris = ±ri2ri4, and consequently the 2kF 
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charge density operator contains either a factor cos(v^(/)o-) (plus sign) or sin(-\/20o-)- A similar 
discrete "gauge covariance" exists of course for all correlation functions. 

In more complicated cases like the two-chain problem, more then one /i-operator occurs. 
Even if they all commute, as is the case for the two-chain problem, additional constraints on the 
permissible eigenvalues of the h's exist due to the existence of relations of the type 

ha/S-yshnX^iuhnpar = il , (A8) 

and similar relations involving more than three h^s. However, a discrete gauge freedom of the 
type mentioned above often remains. For the particular case of fermions with an internal SU(N) 
symmetry, [^] bosonization can be performed without problem and all the /i-operators can be 



given eigenvalue +1, a fact not noticed in the original work. 
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